/  T^J^Va^sSG^v^n&Jf^^^SSVFf JS'WlE^/"^” 


♦■mwhwt **  ’•>«•'**■  *>*  ,‘?n*Jr  ^^Ff^>^^^*^i<^v--; 


AD-752  468 

PROGRAM  SCORES  -  SHIP  STRUCTURAL  RESPONSE 
IN  WAVES 

A  lfred  I .  Raff 

Oceanics,  Incorporated 


Prepared  for: 

Ship  Structure  Committee 
Naval  Ship  Engineering  Center 


July  1972 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


PROGRAM  SCORES— SHIP  STRUCTURAL 
RESPONSE  IN  WAVES 


This  document  has  boon  a&prevod 
for  pybik  roleaso  mi  sate;  Its 
distribution  1$  unlimKod. 


R»*»>r  f . •  i  \  i 

NATIONAL  TECHNICAL 
iNfOR/ZATIGN  SERVICE 


D  D  C 

njETSEHQIiEJ 

'A  DEC  14  1972 


•*  <  v  a  .:;5i 


CGMffim 


t¥n 


KHfP  STRUCTURES  COK38SITTHK 

AH  INTERAGENCY  ADVISORY 
COMMIYTEE  DEDICATED  TO  IMPROVING 
THE  STRUCTURE  OF  SHIPS 


MEMBER  AGENCIES: 


ADDRESS  CORRESPONDENCE  TO: 


UNITED  STATES  COAST  GUARD 
NAVAt  SHIP  SYSTEMS  COMMAND 
MILITARY  SEALIFT  COMMAND 
MARITIME  ADMINISTRATION 
AMERICAN  BUREAU  OF  SHIPPING 


SECRETARY 

SHIP  STRUCTURE  COMMITTEE 
U.S.  COAST  GUARD  HEADQUARTERS 
WASHINGTON.  O.C.  20591 


SR- 174 
1972 


Dear  Sir: 

A  major  portion  of  tl'-e  effort  of  the  Ship  Structure  Committee 
program  has  been  devoted  to  improving  capability  of  predicting 
the  loads  which  a  ship’s  hull  experiences . 

This  report  contains  details  of  a  coaputer  program,  SCORES, 
which  predicts  these  loads.  Details  of  the  development  and 
verification  of  the  program  are  contained  in  SSC-229,  Evaluation 
and  Verification  of  Computer  Calculations  of  Wave-Induced  Ship 
Structural  Loads.  Additional  information  on  this  program  may 
be  found  in  SSC-231,  Further  Studies  of  Computer  Simulation  of 
Slasadng  and  Other  Wave-Induced  Vibratory  Structural  Loadings. 

Comments  on  this  report  would  be  welcomed. 


Sincerely, 


i 


_ unt\  Classification 

'  DOCUMENT  CONTROL  DATA  ■  R  4  D 

J*«  i  «  lattifii  »*l  fitly  6, "Jy  *  I  <ib  fr<n  f  miJ  ntjvutitfi  .utnutnti  n  n  u-»f  l> <  vfifrfd  tthin  fli*  overall  report  i*  < 

1  Qt,  Ji»,»  •  i  N  G  *  C  '  1  V1  f  ^  t  Corporate  author »  RCTOR'  SEC^RUv  CL*s$i»  'CA’iON 

Unclassified 

Oceanics,  Inc.  - 

Plainview,  New  York 

)  «».  P  OR  t  t  .  *  u  F 


Program  SCORES  —  Ship  Structural  Response  in  Waves 


4  ycscniptive  NOTES  (Type  ot  report  and  indtt  u\  e  date  *• » 


[S  a  „  'mQRiJi  fFtrat  nan-e,  n.l'tJI*  initial,  la»t  'mine* 


A.  I.  Raff 


it  REPOR T  CATC 


July  1972 

M  CONf  N*C  ’  OR  Cfi*NT  NO 

N00024-70-C-5076 


6  PROJECT  no 


MC  OTR'fiuTiO^  jrAUwENT 


f  TOTAL  no  Of  paces  76  NO  O  f  REF1 

53 

7a  OPiC'NATOR'S  REPORT  MjMBP  Rl)l 


7b  other  REPORT  NOiS)  (Any  other  number*  that  may  be  aa  tlgned 
thi *  report) 


SSC-230 


Unlimited 


n  supplementary  ,otes 


[d  abstract 


'i  SPONSORING  MLI'AR.  ACTIVITY 


Naval  Ship  Systems  Command 


Information  necessary  for  the  use  of  the  SCORES  digital  compu¬ 
ter  program  is  given.  This  program  calculates  both  the  vertical  and 
lateral  plane  motions  and  applied  loads  of  a  ship  in  waves.  Strip 
theory  is  used  and  each  ship  hull  cross-section  is  assumed  to  be  of 
Lewis  form  for  the  purpose  of  calculating  hydrodynamic  forces.  The 
ship  can  be  at  any  heading,  relative  to  the  wave  direction.  Both 
regular  and  irregular  wave  results  can  be  obtained,  including  short 
crested  seas  (directional  wave  spectrum).  All  three  primary  ship 
hull  loadings  are  computed,  i.e.  vertical  bending,  lateral  bending 
and  torsional  moments. 

All  the  basic  equations  used  in  the  analysis  are  given,  as 
well  as  a  description  of  the  overall  program  structure.  The  input 
data  requirements  and  format  are  specified.  Sample  input  and  out¬ 
put  are  shown.  The  Appendices  include  a  description  of  the  FORTRAN 
program  organization,  together  with  flowcharts  and  a  complete  cross- 
referenced  listing  of  the  source  language. 
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ABSTRACT 


Information  necessary  for  the  use  of  the  SCORES  digital  compu¬ 
ter  program  is  given.  This  program  calculates  both  the  vertical  and 
lateral  plane  motions  and  applied  loads  of  a  ship  in  waves.  Strip 
theory  is  used  and  each  ship  hull  cross-section  is  assumed  to  be  of 
Lewis  form  for  the  purpose  of  calculating  hydrodynamic  forces.  The 
ship  can  be  at  any  heading,  relative  to  the  wave  direction.  Both 
regular  and  irregular  wave  results  can  be  obtained,  including  short 
crested  seas  (directional  wave  spectrum).  All  three  primary  ship 
hull  loadings  are  computed ,  i.e.  vertical  bending,  lateral  bending 
and  torsional  moments. 

All  the  basic  equations  used  in  the  analysis  are  given,  as 
well  as  a  description  of  the  overall  program  structure.  The  input 
data  requirements  and  format  are  specified.  Sample  input  and  out¬ 
put  are  shown.  The  Appendices  include  a  description  of  the  FORTRAN 
program  organization,  together  with  flowcharts  and  a  complete  cross- 
referenced  listing  of  the  source  language. 
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I .  INTRODUCTION 


This  manual  describes  in  detail  the  use  of  SCORES, 
which  is  a  digital  computer  program  for  the  calculation  of  the 
wave-induced  motions  and  loads  of  a  ship.  Both  the  vertical  and 
lateral  plane  motions  are  treated,  so  that  results  for  vertical 
bending,  lateral  bending  and  torsional  hull  moments  can  be  ob¬ 
tained.  The  principal  assumptions  of  the  method  are  that  the 
motions  are  linear,  can  be  solved  by  "strip  theory"  and  that 
the  ship  sections  can  be  approximated  by* "Lewis  forms"  for  the 
purpose  of  calculating  the  hydrodynamic  forces,  that  is,  the 
required  two-dimensional  added  mass  and  wave  damping  properties 
Both  regular  or  irregular  waves  can  be  specified,  and  for  the 
latter  multi-directional  (short  crested)  seas  are  allowed. 

oOORES  was  written  in  thi  FORTRAN  IV  language  and 
checked  out  and  run  on  the  Control  Data  6600  Computer  using  the 
SCOPE  operating  system  (version  3.1.6).  The  program  is  un¬ 
classified. 

The  method  of  analysis  used  in  SCORES  is  outlined  below 
in  Section  II.  All  the  equations  of  motion  and  loadings  are 
given.  In  Section  III,  the  organization  of  the  SCORES  program 
is  discussed  briefly.  An  explanation  of  input  data  card  prepara¬ 
tion  is  given  in  Section  IV,  and  of  program  output  in  Section  V.' 

An  example  problem  is  shown.  Error  messages  which  can  appear 
during  program  execution  are  described  in  Section  VI. 

The  Appendices  include  a  description  of  the  FORTRAN 
program  organization,  flowcharts  for  each  subprogram  and  a  com¬ 
plete  crosp-referenced  (to  the  flowcharts)  listing  of  the  source 
language . 

II.  METHOD  OF  ANALYSIS 

The  analysis  used  in  SCORES  was  developed  and  investigated 
to  some  extent  in  work  supported  by  the  Ship  Structure  Committee.* 

The  exposition  to  be  given  here  will  serve  as  a  convenient  listing 
of  the  equations,  but  for  the  full  derivation  and  explanation  of 
the  analysis  method,  the  references  listed  should  be  consulted. 

•Kaplan,  Paul,  "Development  of  Mathematical  Models  for  Describing 
Ship  Structural  Response  in  Waves,"  Ship  Structure 
Committee  Report  SSC-193,  January  1969  (AD  682591) 

Kaplan,  P. ,  Sargent,  T.P.  and  Raff,A.I.,  "An  Investigation  of  the 
Utility  of  Computer  Simulation  to  Predict  Ship 
Structural  Response  in  Waves,"  Ship  Structure 
Committee  Report  SSC-197,  June  1969  (AD  690229) 

Kaplan,  P.,  and  Raff,  A. I.,  "Evaluation  and  Verification  of  Computer 

Calculations  of  Wave-Induced  Ship  Structural  Response," 
Ship  Structure  Committee  Report  SSC-229,  July  1972, 
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The  relationship  between  the  water  wave  system  and  the 
ship  coordinate  axes  system  is  shown  in  Figure  1.  The  wave  propa¬ 
gation,  at  speed  c,  is  considered  fixed  in  space.  The  ship  then 
travels,  at  speed  V,  at  some  angle,  8  with  respect  to  the  wave 
direction.  The  wave  velocity  potential,  for  simple  deep-water 
waves,  is  then  defined  by: 

<i  =  -ace  cos  k  (x'  +  ct) 
w  (1) 


where  a  =  wave  amplitude 
c  =  wave  speed 
k  =  wave  number  =  — ^ 

A 

\  =  wave  length 

z'  =  vertical  coordinate,  from  undisturbed  water  surface 
positive  downwards 

x'  =  axis  fixed  in  space 

t  =  time 


The  x-y  axes,  with  origin  at  G,  the  center  of  gravity  of  the  ship, 
translate  with  the  ship.  The  x'  coordinate  of  a  point  in  the  x-y 
plane  can  be  defined  by: 

x*  =  -(x+Vt)  cos  e  +y  sin  8  (2) 

Then,  the  surface  wave  elevation  n  (positive  upwards)  can  be  ex¬ 
pressed  as  follows: 


n 


1 

g 


=  a  sin  k  (x*  +  ct) 
z'=o 


(3) 


since 


where  g  =  acceleration  of  gravity 
In  x-y  coordinates  we  have: 


n  =  a  sin  k  [-x  cos  3  +  y  sin  B+(c-V  cos  B)t]  (4) 


•  Dn  ,  3 

n  Dt  3t 


-V  £-)  n  (x,t) 


n  =  akc  cos  k  [-x  cos  3+y  sin  8  +  (c-V  cos  B)t]  (5) 


\ 


I 


II  =  “  6t  [A33  (2“x9+V9)  ]-N^(z>xe+V0)-pgB*(z-X0)  (11) 

where 

p  =  density  of  water 

A^-  local  sectional  vertical  added  mass 

N'  =  local  sectional  vertical  da  oing  force 
coefficient 

E*  =  local  waterline  beam 

and 

N;  =  pg2A^'3  (12) 

with 

A  =  ratio  of  generated  wave  to  neave  amplitude  for 
vertical  motion-induced  wave 


Expanding  the  derivative,  we  obtain: 


dA' 

i  *  -  a33  <*-’«+2ve>  -[  *)i-v^ii 


(z-xe+VG) 


-  ogB* (z-xG)  I 

The  equations  of  motion,  (9)  and  (10)  are  then  transformed  into 
the  familiar  form  as  follows: 


a'z  +  bz  +  c'z  -  de  -  e6  -  g'e  =  Z 


A9  +  B0  +  CG  -  Dz  -  Ez  -  G'z  =  M. 


The  coefficients  on  the  left  hand  sides  are  defined  by: 


a’  =  m+  Al,dx 


b  =  i  N*  dx  -V  d  (Ai,) 
|  2  JJ 


c’=  pg  B*dx 


d  =  D  =  x  dx 


e  =  N'xdx  -2V  A'  dx-V  xd  (A'  ) 
Z  j  J  JJ 


g’  =  pg  I  B*xdx  -Vb 


A  =  1^+  A33  x2dx 


was* 
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B  =  N^x2dx  -2V  A^xdx  -V  x2d  (A^) 


C  =  pg  B*x2dx-VE 


E  =  N^xdx-V  xd  (A^) 


G'=  pg  B*xdx 
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where  ii  =  mean  section  draft.  Substituting  the  expressions  for  n, 
n  and  n  from  Eq.  (4),  (5)  and  (6),  with  y=0  and  applying  the 
approximate  factor  for  short  wave  lengths  we  ohtain: 


(pgB*=A23  kg)sin(-kx  cos  6)  + 


kc  (N^-V  — cos  (-kx  cos  B>  cos  icgt  +  (ogB*-A^kg) 


I  ^33 

cos(-kx  cos  B)-kc  N^-V  “33T”  cos  8)  sin  u  t  J- • 


(20) 


(21) 


where  H  =  local  section  draft 

C  =  local  section  area  coefficient 
s 

The  steady  state  solution  of  the  equations  of  motion  are 
obtained  by  conventional  methods  for  second  order  ordinary 
differential  equations,  using  complex  notation.  The  solutions  are 
expressed  as: 


z  =  z„  sin  (w  t+5) 
o  e 

(22) 

6=6  sin  (u  t+e) 
o  e 

where  the  zero  subscripted  quantities  are  the  amplitudes  and  6 
e  are  the  phase  angle  differences,  i.e.  leads  with  respect  to  the 
wave  elevation  in  Eq.  (7) . 

The  local  vertical  loading  is  given  by: 


azw 

~3x 


(23) 


^M-fl*i,'''«rcitfv^y>v«i^»'  f-  ^**ov  •* 


where 


6m  =  local  mass,  per  unit  length. 


Bq.  (23)  is  simply  the  summation  of  inertial,  hydrodynamic,  hydro* 
static  and  wave  excitation  forces.  The  latter  terms  are  given  in 
Eas.  (13)and  (20).  The  vertical  bending  moment  at  location  xo  is 
then  given  by: 


BMZ  <xo>  «  or 


(x-xQ)  zr  d*  (24) 


and  is  expressed  in  a  form  similar  to  the  motions,  i.e. 


BM„  =  BM  sin  (w  t+o) 
z  zo  e 


B.  Lateral  Plane  Equa.  j.ons 

The  coupled  equations  of  motion  for  sway,  y  (positive  to 
starboard),  yaw,  4  (positive  bow-starboard),  and  roll,  4  (positive 
starboard-down),  are  given  as: 


—  dx+Y 
dx  w 


I  4  “  I  <| 
z  xz1 


§  *  a*+sw 


I  4  —I  4  = 

X*  XZ^ 


||  dx-mg  GM  4+Kw 


where  Iz  *  mass  moment  of  inertia  of  ship  about  z  axis 

Ix  =  mass  moment  of  inertia  of  ship  about  x  axis 

I  =  mass  product  of  inertia  of  ship  in  x-z  plane 
xz 
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Jy 

*  local  sectional  lateral  hydrodynamic  force  on  ship 
ciX 

■=  local  sectional  hydrodynamic  rolling  moment  on  ship 
Yw'  V  Kw  55  wave  excitation  force  and  moments  on  ship 
GM  =  initial  metacentric  height  of  ship  (hydrostatic) . 
The  hydrodynamic  force  and  moment  are  taken  to  be: 


where  OG  =  distance  of  ship  C.G.  from  waterline,  positive  up 
M  =  sectional  lateral  ai  _d  mass 

S 

Ng  =  sectional  lateral  damping  force  coefficient 

M  .  »  sectional  added  mass  moment  of  inertia  due  to  lateral 
89  motion 

Ng.  =  sectional  damping  moment  coefficient  due  to  lateral 
9  motion 

Ir  =  sectional  added  mass  moment  of  inertia 
Nr  =  sectional  damping  moment  coefficient 


Frs  =  sectional  lateral  added  mass  due  to  roll  motion 

Nrs  =  sectional  lateral  damping  force  coefficient  due  to 
roll  motion 

and  the  sectional  added  mass  moments  and  damping  moment  coefficients 
are  taken  with  respect  to  an  axis  at  the  waterline. 


10 


The  additional  roll  damping  moment  to  account  for  viscous 
and  bilge  keel  effects  is  taken  as  a  particular  fraction  of  the 
critical  roll  damping,  as  follows: 

N  *  =  c  C  /L-N  (w  )  (31) 

r  $  c  r  p 

where  Nr*  =  sectional  damping  moment  coefficient  due  to  viscous 
and  bilge  keel  effects 

C  =  fraction  of  critical  roll  damping  (empirical  data) 

<j> 

Cc  =  critical  roll  damping 

L  =  ship  length  (L=x^-Xs ) 

=  natural  roll  (resonant)  frequency 

Nr(w^)  =  va^-ue  Nr  frequency  w  . 

The  critical  roll  damping  is  expressed  in  terms  of  the  natural 
roll  frequency  by: 


where  the  integral  is  over  the  ship  length.  The  calculation  of 
the  natural  roll  frequency,  w  ,  as  indicated  above  is  carried 
out  by  means  of  successive  approximation. 

Expanding  the  derivatives ,  we  obtain 


....  .  |3gg 


The  local  sectional  lateral  force  and  rotational  moment 
due  to  the  waves  acting  on  the  ship  are  represented  as: 


t  f  Dv  dM  Dv  dM  1 

r  -  <»stV  est  -vvw  sr  «8Vk  -Ms«Dt- +v  -as1  vw, 


wB* 

sin  — r—  sm  6 


sin  6 


Dt*Ms*vw* 


B*3  ** 

IT"  "bz 


'w 

— -  -N  v 
t  s$  w 


•  /ttB*  •  0 

sm  sm  8 


sin  8 


where  v  =  lateral  orbital  wave  velocity 
w 

S  =  local  section  area 

2  =  local  sectional  center  of  buoyancy,  from 
waterline 

The  lateral  wave  orbital  velocity  is  obtained  as  follows: 


vw  "  3y 


vw  =  -  akc  e-*^  sin  8  sink  ^-x  cos  8  +  y  sin  6  +  (c-V  cos  8)  tj  (44) 


and  then  we  have: 


w  _  ,  „  -kh  . 

- —  =  -  akg  e  sm 


8  cos  k  jj-x  c.  i  6  +  y  sin  6  +  (c-V  cos  6)tJ  (45) 
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(49) 
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=  BMyo  sin  (wet  +  t) 

TM  =  TM  sin  (u  t  +  v) 
x  xo  e 


(55) 
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The  requirement  on  the  local  vertical  mass  center  is: 

r*b 

6m. ?dx  =  0  (56) 

J*s 

Similarly,  the  requirement  on  the  local  roll  gyradius  is: 

fX. 


6my2dx  =  I 


(57) 


x_ 


The  product  of  inertia  in  the  x-z  plane  is  defined  by: 

rx. 


xz 


6mx!;dx 


(58) 


C.  Wave  Spectra  Equations 

The  wave  spectrum  for  calculations  in  irregular  seas  is 
considered  to  be  a  separable  function  of  wave  frequency  and 
direction  as  follows: 


S  (u,u)  =  (w)  S2(u) 


for  0<w<“ 


(59) 


and  - 


TT  IT 

77  <  U  <  — 
2  —  -  2 


where  S  (co,u)  =  directional  spectrum  of  the  seaway  (short 
crested  sea  spectrum) 


0)  = 


V  = 


S1(io) 

S2(U) 


circular  wave  frequency 

wave  direction  relative  to  predominent  direction 
frequency  spectrum  (long  crested  sea  spectrum) 

spreading  function 


The  SCORES  program  includes  various  spectra  that  can  be 
chosen  as  desired.  However,  in  all  cases,  the  following 
relationship  between  the  spectrum,  or  spectral  density,  and  the 
wave  elevations,  or  amplitudes,  is  used: 


g^TT-v 
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a.,,  =  significant  (average  of  1/3  highest) 

'  wave  amplitude 

ajy^=  average  of  1/10  highest  wave  amplitude. 


Neumann  Spectrum  (1953) 

This  frequency  spectrum  (as  used)  is  given  by: 
S1(«o)  =  0.000827  g21T3u^6e''2g2‘‘,  2u  2 


where  U  =  wind  speed 

The  constant  is  one  half  that  originally  specified  by 
Neumann  so  that  this  spectrum  satisfies  Eq.  (62).  Thus,  originally 
the  Neumann  spectrum  required  only  a  factor  of  /T-  in  Eq.  (65), 
instead  of  2.0.  : 


Pierson-Moskowitz  (1964) 
This  is  given  by: 


S^w)  =  0.0081  g2o)"5e“*74g4a>  4°  4 
and  was  derived  on  the  basis  of  fully  arisen  seas. 


Two  Parameter  (1967) 


S^w)  =  A*  B(u_5e-— u  4 


where 


A  =  0.25  H. 


B  =  (0.817  |^>4 


1^/2  =  significant  wave  height  (=2.0ajy.j) 


T  =  mean  wave  period 


This  spectrum  is  usually  used  in  conjunction  with  "ooserved" 
wave  height  and  period,  which  are  then  taken  to  be  the  significant 
height  and  mean  period.  This  spectrum  is  similar  to  that  adopted 
by  the  I.S.S.C.  (1967)  as  "nominal",  except  that  it  is  expressed 
in  circular  wave  frequency  instead  of  frequency  in  cycles  per 
second. 


*********** .  *  -  I  w»r»^r-rq^>aAij^^ 


?WP^ I^^^f^?5^ES5^!lB[J®55» !  %?AW  7  \K-  "Y~-  'W* 


J&iHSi  fx*r»sn*f*it0itran- a,*^  x 


19 


Uni -Directional  Spreading  (Long  Crested  Seas) 
This  is  obviously: 


Sjfu)  =  3 (p)  (delta  function) 


(70) 


Cosine-Squared  Spreading 


S2(„)  - 


2  ? 
—  cos'p 


(71) 


Responses 


All  of  the  motions  and  moments  calculated  are  considered 
to  be  linear  and  the  principle  of  wave  superposition  is  assumed. 
Thus  for  each  response  a  spectrum  is  calculated  by: 


3^  (w,p)  =  [*L  (u,  p)J  S  (w,  u) 


(72) 


where  T^(w,p)  =  response  amplitude  operator  (amplitude  of  response 


per  unit  wave  amplitude) 

We  then  have,  similar  to  the  wave  amplitude: 


ai 


IT 

\2 


(u,u)  do  dp 


"  | 

8 

1 

J 

rJ 

_ ! 

= 

s2M 

(w,p)J  S]_  (w)  dw 

, 

7T 

*x 

J 

o 

dp 


(73) 


where  a^2  =  mean  squared  response  amplitude. 


Eqs.  (63)  -  (66)  then  apply  to  each  response. 
D.  Non-dimensional  Forms 


Frequency  parameter:  =  — — 

z  g 


H 


1 


Jyf^j 


'i>,  is  ^J^wf  - 

?3BSy 5!Hf^flSj ¥►&*  *  r*  •<  ■>  wot  tw^, 
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Non-dimensional  linear  motion  (heave,  sway):  option  e 


Non-dimensional  angular  motion 
(pitch,  yaw,  roll): 


motion  ampJitude 
2na/X 


Ncn-dimensional  moment: 


BM  (orBM  or  TM  ) 
z  y  x 

og  B*l7a 


Non-dimensional  shear:  . A°Lce 

pg  B*La 

III.  PROGRAM  ORGANIZATION 

A.  General 

In  general,  the  SCORES  computer  program  has  been  arranged 
and  organized  to  both  keep  a)  the  coding  simple  and  flexible  (for 
possible  future  modification)  and  b)  the  running  times  low  (for 
obvious  reasons).  Thus,  precision  of  computation  has  not  been  of 
major  priority  in  program  development.  This  approach  is  considered 
reasonable  at  the  present  time  because  precise  correlation  (to 
less  than  about  5%)  with  independent  data  (model  or  full-scale  ex¬ 
periments)  is  not  envisioned,  and  the  theoretical  analysis  itself 
is  an  approximation. 

Aside  from  the  actual  coding  and  data  structure  in  the 
program,  which  will  not  be  discussed  here  (see  Appendices  A,  B 
and  C  of  this  report) ,  this  approach  manifests  itself  primarily 
in  two  aspects.  The  first  is  the  precision  with  which  the  local,  or 
two-dimensional,  sectional  adj.ed  mass  and  damping  characteristics 
or  properties,  are  calculn.^d.  For  vertical  oscillation,  the  method 
of  Grim*  is  used.  For  the  two-dimensional  properties  in  lateral 
and  roll  oscillations,  the  method  of  Tasai**  has  been  programmed. 

In  general,  these  methods  can  be  carried  out  to  increasing  degrees 
of  numerical  accuracy.  For  practical  purposes  of  keeping  running 
time  reasonable,  these  calculations  have  been  limited.  For  example 
in  the  lateral  and  roll  computations,  the  infinite  series  of  terms 
representing  the  velocity  potential  is  truncated  to  nine  terms 
and  only  15  points  along  the  Lewis  form  contour  are  used  for  least 
square  approximation  purposes.  While  the  full  range  of  section 
properties  and  frequencies  has  not  been  explored  in  detail,  results 
on  the  order  of  1%  accuracy  or  better  are  obtained  for  average 
sections  over  a  wide  frequency  range. 

*  Grim,  0.,  "Die  Schwingungen  von  schwimmeden,  zweidimensionalen 
Korpern,"  Ilf, <?h  Report  No.  1171,  September  1959. 

Grim,  O.,  and  Kirsch,  M. ,  private  communication,  September  1967. 

**Tasar,  F.,  "Hydrodynamic  Force  and  Moment  Produced  by  Swaying  and 
Rolling  Oscillation  of  Cylinders  on  the  Free  Surface," 
Reports  of  Research  Institute  for  Applied  Mechanics, 
Kyushu  University  Japan,  Vol.  IX,  No.  35,  1961 
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The  second  aspect  of  program  organization  is  related  to  the 
above.  While  the  computations  ot  the  two-dimensional  properties 
are  limited  as  described,  they  still  are  relatively  lengthy.  That 
is  at  a  particular  condition  of  ship  speed,  wave  angle  and  wave 
length,  the  bulk  of  the  computation  time  would  be  devoted  to  these 
calculations  rather  than  the  formation  of  the  coefficients, 
wave  excitation,  solution  of  ship  motions  and  the  resulting 
calculation  of  applied  moments.  Therefore,  it  was  decided  that 
rather  than  calculate  for  each  frequency  at  each  cross-section 
the  above  mentioned  two'-dimensional  properties,  instead  the  two- 
dimensional  properties  are  calculated  first  at  25  values  of 
frequency  over  a  wide  range  and  then  interpolated  (or  extra¬ 
polated)  for  each  subsequent  frequency.  The  results  of  the  initial 
calculation  over  the  frequency  range  are  saved  in  the  computer 
memory  for  the  calculations  at  hand,  and  can  also  be  saved  on  a 
permanent  disc  file  (or  magnetic  tape  storage) ,  for  later  usage. 

In  chis  way,  a  large  range  of  ship  speeds  and  headings  car.  be  run, 
each  over  the  appropriate  frequency  range,  without  excessively1 
high  running  times.  The  interpolation  procedure  used  is  a 
six-point  continued  fraction  method1  which  giv^s  results  that  are 
generally  well  within  1%. 

In  other  respects,  the  SCORES  program  is  organized  in  a 
fairly  straightforward  manner.  The  input  consists  of: 

a)  basic  data  which  specify  the  hull  form  and  weight 
distribution  and 

b)  conditional  data  which  specify  the  speeds  and  wave 
parameters . 

Repeated  sets  of  conditional  data  can  be  run  with  the  same  basic 
data,  that  is,  for  the  same  defined  ship.  A  fair  amount  'of  input 
data  verification  is  incorporated  into  the  program. 

B.  Restrictions 

The  main  restrictions  in  the  program  concern  the  following 

items: 

Maximum  no.  of  ship  cross-sections . 21 

(stations  0  to  20) 

Maximum  no.  of  wave  angles  (in  one  run)...,.  .25 

Maximum  no.  of  wave  lengths  (in  one  run).... 51 

Maximum  no.  of  sea  states  (in  one  run) . 10 

The  core  storage  requirement  is  about  25,000  cells  as 
compiled  on  the  CDC  6600.  This  includes  the  program  instructions, 
data  storage  and  system  routines  to  handle  input-output  system 
control  and  provide  mathematical  functions.  It  would  be  possible 
to  decrease  this  core  requirement  via  program  overlay  and 
linkage  techniques,  should  the  need  arise.  However,  it  probably 
would  be  relatively  difficult  to  fit  the  program  within  a  12K 
core  restraint. 


-v-  av  - 
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The  word  length  on  the  CDC  6600  is  60  bits.  No  loss  in 
overall  computational  accuracy  would  be  expected  if  this  were 
reduced,  as  in  other  digital  computers,  to  36  bits. 

A  special  system  subroutine  called  DATE  is  used  which 
provides  the  current  -“ate.  This  is  used  only  in  the  heading  on 
the  output. 

C.  Running  Time 

The  following  approximate  times  are  for  running  under  the 
SCOPE  operating  system  on  the  CDC  6600  computer. 

Program  compilation  (RUN  compiler) . 10. 0  secs. 

Program  loading  into  core .  1.0  secs. 


Calculation  of  TDP*  Array  (21  sections, 
both  vertical  and  lateral  modes) . 


25  secs. 


Calculate  motions,  moments  at  one  condition, 

(21  sections,  both  vertical  and  lateral 

modes) . . . . . . . 0.14  secs. 

Calculate  spectral  response,  for  each 

spectrum,  for  each  condition . .  0.006  secs. 

Thus,  for  a  run  with  two  ship  speeds,  7  headings  (at  30°  increments 
from  head  to  following  seas),  21  wave  frequencies  (to  adequately 
cover  the  spectral  energy  bands)  and  5  sea  states,  the  incremental 
time  once  the  program  was  compiled,  loaded  and  the  TDP  Array  was 
calculated,  would  be  estimated  as  follows: 


(21)  [0 . 14+  (5)  (0.006)]  =  50  secs. 


IV .  DATA  INPUT 


This  section  of  the  manual  describes  the  details  of  data 
card  input  to  the  SCORES  program. 

A.  Units 


For  calculations  in  regular  waves,  there  are  no  inherent 
units  assigned  to  any  of  the  variables  in  the  program.  Thus,  the 
user  is  free  to  choose  any  desired  set  as  long  as  they  are 
consistent  for  all  input  parameters.  The  units  are  established 
by  the  input  values  of  water  density  and  gravity  acceleration. 
Some  typical  units  are  shown  below. 


♦Two-dimensional  properties 
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Water  Density 


Gravity  Accel. 

Resultant  Unit 
System 


lbs./cu.  ft. 


ft. /sec. 2 
ft. -lbs. -sec. 


tons/cu.  ft. 


ft. /sec. ' 


ft . -tons-sec. 


metric  ton/cu. 
meter _ 

meter/sec. 2 

meter-metric 

ton-sec. 


Wave  direction  angles  are  always  specified  in  degrees, 
rather  than  radians . 

However,  for  spectral  calculations  in  irregular  waves,  using 
either  the  Neumann  or  Pierson-Koskowitz  spectra,  the  SCORES  pro¬ 
gram  assumes  ft. -sec.  units,  full  scale.  The  input  wind  speeds 
used  to  specif}'  spectral  intensities,  or  sea  states,  are  then 
assumed  to  be  in  knots. 

The  following  input  data  description  indicates  typical 
consistent  units  for  all  parameters.  Other  systems  of  units 
could  be  used,  as  noted  above. 

B .  Data  Card  Preparation 

Every  data  card  defines  several  parameters  which  are 
required  by  the  program;  each  of.  these  parameters  must  be  input 
according  to  a  specific  format.  "I"  format  (integer)  means  that 
the  value  is  to  be  input  without  a  decimal  point  and  packed  to 
the  right  of  the  specified  field.  "F"  format  (floating  point) 
requires  that  the  data  be  input  with  a  decimal  point;  the  number 
can  appear  anywhere  in  the  field  indicated.  "A"  format 
(alphanumeric)  indicates  that  certain  alphabetic  characters  or 
title  information  must  be  entered  in  the  appropriate  card  columns . 

If  the  field  is  left  blank  for  either  "1"  or  "F"  format, 
a  value  of  zero  (0)  is  assigned  to  the  parameter.  Thus,  parameters 
not  required  by  the  program  for  a  particular  problem  need  not  be 
specified. 

The  card  order  of  the  data  deck  must  follow  the  order  in 
which  they  are  described  below.  Cards  which  must  be  present  in 
every  run,  regardless  of  options,  are  marked  with  an  asterisk  (*) . 
The  first  eight  types  of  cards  are  considered  the  basic  data  set, 
while  subsequent  cards  are  the.  conditional  data  set(s). 

1)  Title  Card  (*} 


Columns 


Format 


Entry 

Ar.y  alphanumeric  title 
information,  used  to  label 
job  output 


\a&esawmMLH 


The  first  30  columns  are  used  a3  a  label  for  the  TDP  array  file. 
Thus,  subsequent  runs  using  the  file  must  duplicate  these  first 
30  columns  which  are  then  checked  against  the  file  label  before 
using  the  data.  This  avoids  inadvertent  use  of  an  incorrect 
TDP  file. 

2)  Option  Control  Card  (*) 


Columns 

1-2 

3-4 

5-6 

7-8 

9-10 

11-12 

13-14 

15-16 

17-18 

19-20 

21-22 


Format  Entry 

I  Integration  option  control  tag 

I  Moment  option  control  tag 

I  Mass  dist.  option  control  tag 

I  Wave  spectra  option  control  tag 

I  Degrees  of  freedom  option  control  tag 

I  Directionality  option  control  tag 

I  TDP  file  option  control  tag 

I  Moment  closure  option  control  tag 

I  Output  form  option  control  tag 

I  Torsion  axis  option  control  tag 

I  Number  of  ship  segments 


Each  option  control  tag  is  given  a  value  of  0,  1,  2  or  3 
where  the  meaning  of  each  is  given  in  the  table  below.  The  last 
entry  of  the  card,  the  number  of  ship  segments,  corresponds  to  the 
even  number  of  equal  length  segments,  or  strips,  into  which  the 
ship  hull  is  divided  lengthwise  for  purposes  of  calculation. 


OPTION  CONTROL  TAG  INTERPRETATION 


Letter  | 
Code  ■ 


Tag 

Descriptor 


Options  Available 


B  i 


Integration 


Moment 


0:  Simple  summation 
Is  Trapezoidal  rule 

0:  Calc,  motions  only,  use 

summary  mass  properties 
1:  Calc,  motions  only,  use 
mass  dist. 

2:  Calc,  moments,  use  mass 
dist. 


Mass  dist. 


0:  Input  masses 
1:  Input  weights 


Wave  spectra  0:  Regular  waves 

1:  Neumann  spectra 
2:  Pierson-Moskowitz  spectra 
3:  Two  parameter  spectra 

(continued  on  next  page) 


OPTION  CONTROL  TAG  INTERPRETATION,  Continued 


Letter 

Code 

Tag 

Descriptor 

[  Options  Available 

E 

i 

Degrees  of 
freedom 

i 

0:  Vertical  plane  only 

1:  Vertical  and  lateral  plai 
2:  Lateral  plana  only 

F 

Direction¬ 

ality 

0:  Uni-directional  waves 

Is  Cos-sq.  wave  spreading 

G 

1 

TDP  file 

i 

0:  Generate  TDP  file,  write 
on  file  (Tape  10) 

Is  Read  TDP  file, (Tape  10) , 
out  TDP  data 

2:  Read  TDP  file, (Tape  10), 
print-out 

H 

Moment 

closure 

0:  Suppress  closure  calcs. 

Is  Calc,  and  print  out 
closure  results 

I 

Output  form 

0s  Dimensional 

Is  Non-dimensional 

J 

Torsion  axis 

0s  Center  of  gravity 

Is  Waterline 

3)  Length  Card  ( * ) 


Columns  Format  Entry 


11-20  F 
21-30  F 
31-40  F 
41-50  F 


Ship  length  (ft.) 

Water  density  (tons/cu.ft. ) 
Acceleration  of  gravity  (ft. /sec.2) 
Ship  displacement  (tons) 


The  entries  on  this  card  are  self  descriptive  and  determine 
the  units  to  be  used  for  all  other  parameters ,  except  as  noted 
earlier. 


4)  Hull  Form  Cards  (*) 


Columns  Format  Entry 


1-10  F 

11-20  F 

21-30  F 

31-40  F 


Section  waterline  breadth  (ft.) 
Section  area  coefficient  (-) 
Section  draft  (ft.) 

Section  centroid  (ft.) 


u  1.^54* 
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One  card  is  used  for  each  section  to  be  specified,  in  order 
along  the  ship  length  starting  at  the  bow.  For  example,  if  the 
number  of  segments  is  10,  and  the  integration  option  tag  is  0, 
then  10  hull  form  cards  are  required  which  correspond  to  the  hull 
at  stations  1/2,  1  1/2,  2  1/2,  . ..,  8  1/2,  9  1/2.  If  the 
integration  tag  is  1,  then  11  hull  form  cards  are  required  at 
stations  0,  1,  2,  3  .....  9,  10. 

The  entries  for  sectional  waterline  breadth,  area  coef¬ 
ficient  and  draft  are  straightforward.  The  fourth  entry,  the 
section  centroid,  is  measured  downwards  from  the  waterline  .  If 
no  entries  are  given  and  the  centroids  are  needed  for  lateral 
plane  motions  calculations,  apDroxirate  controids  are  then 
calculated  based  on  the  area  coefficient  and  draft  (using  a  two- 
dimensional  version  of  the  Moorish  Approximation) . 

5)  Lateral  Plane  Card 

Columns  Format  Entry 

1-10  F  Ship  vertical  center  of  gravity  (ft.) 

11-20  F  Radius  of  gyration  in  roll  (fv. .) 

This  card  is  used  only  if  the  degrees  of  freedom  option 
tag  is  1  or  2,  indicating  lateral  plane  calculations.  The  ship 
vertical  c.g.  is  measured  from  the  waterline,  positive  upwards. 

6)  Summary  Mass  Properties  Card 

Columns  Format  Entry 

1-10  F  Radius  of  gyration,  longitudinal 

(ft.) 

11-20  F  Longitudinal  center  of  gravity 

(ft.) 

This  card  is  used  only  if  the  moment  option  tag  is  0. 

The  longitudinal  center  of  gravity  is  measured  from  amidships, 
positive  forwards. 


7)  Sectional 

Mass 

Properties  Cards 

Column 

Format  Entry 

1-10 

F 

Segment  weight,  or  mass  (tons 
or  tons-sec2/ft. ) 

11-20 

F 

Segment  vert.  c.g.  (ft.) 

21-30 

F 

Segment  roll  gyradius  (ft.) 

These  cards  are  used  only  if  the  moment  option  tag  is 
1  or  2,  in  lieu  of  the  summary  mass  properties  card  above.  One 
card  is  used  for  each  section  to  be  specified,  in  a  similar 
manner  as  the  hull  form  cards  described  earlier. 

The  first  entry  on  each  card  is  the  segment  weight,  or 
mass,  depending  on  whether  the  mass  dist.  option  tag  is  1,  or  0, 
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respectively.  The  second  entry,  the  segment  vertical  center  of 
gravity,  necessary  only  for  lateral  bending  moment  calculations, 
is  measured,  positive  downwards,  with  respect  to  the  ship's  over¬ 
all  vertical  center,  as  specified  on  the  lateral  plane  data  card 
above.  Since  it  is  required  that  the  vertical  mass  moment 
integral  satisfy  the  specified  overall  v.c.g.,  the  input  segment 
v.c.g.'s  are  shifted  by  an  equal  amount,  up  or  down  as  necessary 
to  exactly  balance  the  vertical  moment  for  the  hull.  This 
minimizes  the  effort  required  to  obtain  precise  balance  in  input 
data  preparation.  The  third  card  entry,  the  segment  roll  gyradius, 
is  needed  only  for  torsional  moment  calculations.  If  no  entries 
are  given  the  overall  ship  value  is  used  at  each  segment. 

8)  Moment  Station  Card  (*) 

Column  Format  Entry 

1-10  1  First  station  for  moment  calculations 

11-20  I  Last  station  for  moment  calculations 

21-30  1  Increment  between  stations 

The  parameters  on  this  card  determine  where  along  the  ship 
hull  the  moment  calculations  are  to  be  performed.  Station  numbers 
are  defined  as  zero  at  the  forward  end  of  the  first  segment, 
increasing  to  N,  the  number  of  segments,  at  the  after  end  of  the 
last  segment.  If  the  calculations  are  required  only  at  one  station, 
then  the  first  two  entries  on  the  card  snould  be  equal  to  that 
station  number. 

The  moment  results  at  only  one  station  are  stored  for 
subsequent  irregular  seas  spectral  calculations.  In  the  calculations 
over  a  range  of  stations  at  which  moments  are  calculated  (and 
printed) ,  then  only  the  results  at  midships  are  stored  for  the 
subsequent  spectral  calculations. 


9)  Run  Control  Card  (*) 


Columns 


Format 


1-10  F  Run  control  tag  and  wave 

amplitude  (ft.) 

11-20  F  Initial  wave  length,  or 

frequency  (ft.  or  rad. /sec.) 

21-30  F  Final  wave  length,  or  frequency 

(ft.  or  rad. /sec.) 

31-40  F  Increment  in  wave  length,  or 

frequency  (ft.  or  rad. /sec.) 

41-50  F  Initial  ship  speed  (ft. /sec.) 

51-60  F  Final  ship  speed  (ft. /sec.) 

61-70  F  Increment  in  ship  speed  (ft. /sec.) 

The  first  entry,  the  run  control  tag,  determines  program 
continuity: 


«****«*« 


Run  Control  Tac 


Action 


Greater  than  0.0 


Continue  calculations,  using  this  as 
wave  amplitude 


0.0 (or  blan) 


Stop  calculations;  read  new  basic 
data  set 


Less  than  0.0 


Stop  program  execution 


Thus,  if  the  run  control  tag  is  not  greater  than  0.0,  then 
the  remaining  parameters  on  the  card  are  irrelevant.  A  blank 
card,  for  example,  is  used  to  stop  calculations  and  proceed  to 
read  a  complete  new  set  of  data  starting  with  the  title  card  , 

1)  above.  This  parameter  is  also  used  as  the  wave  amplitude,  and 
is  usually  set  to  1.0. 


The  next  three  entries  determine  the  wave  lengths  to  be 
used  in  the  calculations.  If  the  wave  spectra  option  control  tag 
is  0,  indicating  regular  waves,  then  these  entries  are  the  initial, 
final  and  increment  in  wave  length.  If  the  wave  spectra  option 
control  tag  is  greater  than  0,  indicating  irregular  wave  calculations, 
then  these  entries  are  the  initial,  final  and  increment  in  wave 
frequency.  The  increments  should  always  be  positive,  so  that  wave 
length,  or  frequency,  increases  from  initial  to  final  value. 


The  last  three  entries  sre  similar  parameters  for  ship  speed. 
If  calculations  are  required  at  only  one  value,  then  the  initial 
and  final  values  should  both  be  set  equal  to  it. 


Roll  Dc 


Column 


Format 


Fraction  of  critical  roll  damping 
(empirical  data) 


This  card  is  used  only  if  the  degrees  of  freedom  option 
control  tag  is  1  or  2  indicating  lateral  plane  motions  calculations 
are  included.  The  calculated  wave  damping  in  roll,  at  the  natural 
roll  frequency,  is  increased  so  that  the  total  damping  is  the 
specified  fraction  of  critical  damping.  The  additional  roll 
damping  thus  determined  initially  is  then  used  for  all  subsequent 
calculations . 


11)  Wave  Angle  Card  (*) 


Column 


Format 


Entry 


1-10 

11-20 

21-30 


Initial  wave  angle,  degrees 
Final  wave  angle,  degrees 
Increment  in  wave  angle,  degrees 


These  entries  specify  the  wave  direction  angles  to  be  used 
in  the  calculations  and  are  always  given  in  degrees.  For 
calculations  with  uni-directional  waves,  the  meaning  of  the 
parameters  is  as  indicated.  If  the  directionality  option  control 
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tag  is  greater  than  0,  indicating  calculations  for  a  directional 
wave  spectrum,  then  only  two  choices  exist.  If  the  initial  wave 
angle  is  180.0  the  calculations  proceed  for  head  seas  only, 
including  the  wave  directionality.  If  the  initial  wave  angle  is 
not  180.6  the  calculations  proceed  for  all  angles  from  following 
seas  to  head  seas ,  in  s',teps  according  to  the  wave  angle  increment 
specified. 

In  both  cases  the  integrations  with  respect  to  wave  angle 
use  the  same  increment,  as  specified. 


Wave  Spectra  Card(s) 


Columns 

Format 

Entry 

1-10 

I 

No.  of  sea  states  (wave  spectra) 

11-15 

F 

First  spectra  parameter 

16-20 

F 

Second  spectra  parameter 

21-25 
(5  col. 

F 

Third  spectra  parameter 

fields) 

F 

: 

56-60 

F 

Tenth  spectra  parameter 

This  card  is  used  only  for  calculations  in  irregular  seas 
(wave  spectra  option  control  tag  is  greater  than  0) .  The  first 
entry  specifies  the  number  of  sea  states  (spectra)  to  be  used 
(maximum  10) .  For  both  the  Neumann  and  Pierson-Moskowitz  spectra 
(wave  spectra  option  control  tag  equals  1  or  2),  the  parameters 
to  be  specified  are  the  wind  speed,  in  knots,  for  each  sea 
state.  For  the  two  parameter  spectrum  (option  tag  equals  3), 
tue  parameters  on  this  card  are  the  significant  wave  heights 
for  each  sea  state.  A  second  card  is  then  used  which  contains 
the  mean  periods  for  each  corresponding  sea  state,  as  the 
spectral  parameter  entries  specified  above. 


Sample  Input  Deck 


A  sample  input  card  deck  listing  is  given  on  the  next 
page.  The  units  are  meters,  metric  tons  and  seconds. 


V.  PROGRAM  OUTPUT 


Description 


The  printed  output  from  the  SCORES  program  depends  on  the 
option  control  tags  set  as  input.  Each  output  section  will  be 
described,  though  in  any  given  run  not  all  sections  will  be 
printed.  Each  section  starts  a  new  page  and  is  labeled  with  the 
title  information  and  date. 


The  first  part  of  the  output  is  a  listing  of  the  basic 
input  data  as  processed.  This  defines  the  hull  form  and  weight 
distribution.  Then  the  conditional  data  cards  are  printed  out. 
For  irregular  seas  cases,  the  wave  spectra  will  then  be  printed, 
together  with  internally  generated  wave  statistics.  If  the  TDP 
array  is  calculated  diagnostic  messages  concerning  these 
calculations  may  then  appear. 
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The  next  output  will  be  the  listing  of  the  two-dimensional 
properties  (TDP  array)  for  each  station  and  each  frequency.  If 
the  data  is  being  read  from  file,  this  output  can  be  suppressed. 
For  lateral  plane  calculations,  the  natural  roll  frequency  and 
roll  damping  information  will  then  be  printed. 

Then,  the  vertical  and/or  lateral  plane  responses  will  be 
printed  out  with  all  frequencies,  or  wave  lengths,  for  a  given 
ship  speed  and  wave  angle,  on  the  same  page.  For  irregular  seas 
calculations,  this  will  be  followed  by  a  print-out  of  the 
response  spectra  and  statistics  (long  crested  seas) .  These  pages 
will  be  repeated  for  each  wave  angle  at  the  initial  ship  speed. 
Then  directional  seas  calculations  results  will  be  output,  if 
specified.  The  output  is,  of  course,  then  repeated  for 
additionally  specified  ship  speeds. 

B.  Sample  Output 

A  sample  output  listing,  in  abbreviated  form,  is  given 
following  the  sample  input  listing. 


Sample  Input  Card  Deck  Listing 


StHIFS  6< 

HULL 

•  0.  «o 

12  13 

1  0  1  1  l 

l^j.n 

l.W 

oo.oo 

.0 

rO.W 

1*.39 

ll.OJ 

22.88 

,A*4 

11.03 

26.58 

•  «*9 

11.93 

27.5* 

.9/0 

11. 

27.57  " 

.<*V| 

11.93 

27.57 

11.03 

27.57 

,OV4 

11.03 

27.57 

*9V4 

11.03 

27.57 

•  9V4 

11.03 

27.57 

.994 

11.03 

27.57 

,994 

IU03 

27.57 

#9V1 

11.03 

27.57 

,9»9 

11.03 

27.57 

.9t>* 

11.03 

27.2* 

.9*1 

11.03 

25.9* 

11.03 

231*6 

•  Tb«J 

11.03 

19,63 

.0*7 

11.03 

13.87 

.419 

11.03 

*.«1 

.53 

1.19 

-1.0985 

2*0.  6 
'*81.3 
1203.2 
2*06.3 

3050. 1 

♦ooo.r 

*331.* 

*33l  ,* . . 

3368.6 
168*,* 

168*.* 

1**3, 8 
2198.8 
"3290, 7 

3633.6 
3*4b.t 
31*6, J 

1 935. 1 
121.9 

*«T".3 - - 

120.3 

"  10  ~Ifl  " 

t.O  0.3157  1.307V 

o.io 

10.0  170,0  20.0 

■■ "  - 1  o.* . 

10.0 


BLOC*  (7X0  MOT »  VO,  100  8)  OCt'OVICS  I’Mfl  IfCT  90,  1003 
9.80665  *8126. n 


1 

O.OoSl  6.5257  6.9257  1.0 


-1.0 


I 
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Sample  Input  Listing 


SEMES  NO  HULL  FORM,  O.Hfl  PLUCK  UNO  HP!.  HO.  100  SI  OCEANICS  PPOJECT  NO.  1093  SEP  ?*•  1970 

OPTION  CONTNOL  TAGS  -APCOEFGHtJ 
- -  - 1*  ?~1  3  1  ~0  1  1  1  1  *  WO.  or  STATIONS  »  80 


BASIC  INPUT  DATA 


LENGTH  ■ 

193.00 

density  ■  l. 

025000 

OlSf-w*  ■ 

_4bl?6.40 

GRAVITY  a  9, 

806650 

STATION 

Ht  AH 

AREA  C0fc>, 

DRAFT 

1-ftAP 

HEIGHT 

zfta 

SrR.ROLL 

0*00 

0,0000 

0.OOU0 

0,0000 

0C 000ft 

240,6000 

0,0000 

8.9602 

woo 

14.3900 

.8720 

11  .uJOO 

5.044* 

4ttl .3000 

0.0000 

8,9602 

2  •  00 

22 , flhOO 

,  A940 

11.0300 

5. l?Si 

1203.2000 

o.oooo 

8.9602 

3*00 

28.5800 

.9290 

11.0300 

5.254ft 

2406.3000 

0.0000 

8.9602 

~  4.00 

27,5400 

.9700 

11.0300 

5.4049 

'3B5o;ioo<r 

O.DOOO 

■  8.9602* 

5,00 

27.5700 

.9910 

11.0300 

5.481o 

4090,7000 

d.0000 

8.9602 

6,00 

27.5700 

.9940 

11.0300 

5.492o 

4331.4000 

o.oooo 

8.9602 

7,00 

27.5700 

,9940 

11.0300 

S.492o 

4331.4000 

0.0000 

8.9602 

8.00 

27.5700 

•  9940 

11.0300 

5.492o 

3368.8000 

0.0000 

8,9602 

9,00 

27,5700 

.9940 

11.0300 

5.492o 

1684,4000 

0.0000 

8,9602 

“'*10,00 

27.5700 

“*.9940 

11.0300 

5.4926 

1684.4000 

■  0.0000  " 

8.9602  -  - 

U,00 

27.5700 

.9940 

11.0300 

5,492a 

1443.8000 

0.0000 

8.9602 

12.00 

27.5700 

.9910 

11.0300 

5.489i 

2195.8000 

0.0000 

8,9602 

13,00 

27,5700 

.9890 

11.0300 

5.474* 

3290.7000 

0,0000 

8.9602 

14,00 

27.5700 

,9680 

11.0300 

5.397i 

3633,6000 

0,0000 

8.9602 

15.00 

27,2400 

.9210 

11.0300 

5.224k 

3465.1000 

0.0000 

8.9602 

“16,00 

25.9400 

.8610 

11.0300 

"  4,967? 

3146.3000 

■'  0,0000 

8.9602 

17.00 

23.4600 

.7580 

11.0300 

4.625? 

1955.1000 

0,0000 

8,9602 

■  Ift.no 

1 9. *300 

.6270 

11.0300 

4.143* 

721.9000 

0.0000 

8.9602 

19,00 

13,8700 

.4190 

11.0300 

3.378o 

481 .3000 

0.0000 

8.9602 

20.00 

4.4100 

.5300 

I. 1000 

.3777 

120.3000 

0.0000 

8.9602 

OS  •  -1.099  c.vHAOTUS'.ROU  »  B.9S0 

CALCULATE  MOMENTS  AT  STATION  10 


UtMlVtU  RESULTS 

-ISPL.IWTS.I  •  *8126.50 

-  cons.-  'cur."  f  n.nnrwftr  »rosuiPsr  •  nispl.  ivol.)  *  a*ott.s3 

LONS,  C.O.  s  A.«?S  (EMO.  OF  MlOSMIPSI  LnNO,  OYPAOIUS  *  Aft. 159  OM  n  1.178 


SERIES  NO  •  uL  FOHm,  0.80  HLOCK  (TnO  PPI,  WO.  100  SI  OCEAWICS  PPOJECT  NO.  1093  SFP  ?*•  1970 


CONDITIONAL  INPUT  nATA  CARD  PRINT  OUT 

1.0000  .3157  1,3079  .0A91  6.c257  6,9257  1.0000 

.1000 

_  10,0000  170.0000  80,0000 

1  «,»  -0.0  -0.0  -0.0  -0,0  -0.0  -0,0  .0,0  -0.0  -0.0 

. . ”10.0  -0.0  -o.o-  PP.O  -0.0  -0.0  -0.0  .0.0  -0.0  -0.0 


SERIES  NO  MULL  FORM,  0.H0  BLOCK  UNO  RPT .  NO.  100  SI  OCEAWICS  PROJECT  NO.  1093  SEP  ?*•  19T0 

WAVE  SPECTRAL JiENSJIV*  Tko  PARAMETER.  ISSC  1967  SPECTRA  _  _  _ _ _ 

SIfi.HT.  P.AOO 

Mw, PER .  In, 000 

SPECTRA  NO.  1 

MAVt  FREQ. 

731V  . ."W 

.361  3.328 

.AON  B.NlO 

♦A51  12.254 

1 *A9N  12.954 

.541  11.7*3 

"  “  .9*6 - 9.624 

.631  7.886 

.676  N.20N 

.72?  4. 846 

.767  3, TO? 

.812  ?,96! 

‘  .R5f  '  ?.331 

.902  1.8*7 

.9*7  1 ,*75 

.992  1.18N 

1.037  .961 

1.082  .78* 

1.127  '  .6*4' 


1.173  .533 
1.218  .**3 
1.263  ,371 
1.308  ,313 

- KW7-SU - “A-.ZW 

R.M.S.  2.073 
AV8.  2,539 
SJG.  4.1*6 
AVl/10  5,277 


-  "-'■gns^gjjwis^v  'JS* '  " '  £&*~y**F&5sZ**&7? ?fr37‘- 


^rj^5p^ar» 
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VI.  ERROR  MESSAGES 


Various  error  messages  may  appear  in  the  output  and  cause 
program  termination.  Each  will  be  labeled  with  the  subroutine 
which  found  the  error,  and  possibly  a  brief  note  as  to  the  type 
of  error.  Some  messages  give  error  numbers  as  explained  below: 


Subroutine 

Error  No. 

Explanation 

PRELIMB/C 

0 

Too  many  sections,  wave 
lengths,  wave  angles,  etc. 

PRELIMB 

1 

Sum  of  weight  distribution^ 
displacement 

PRELIMB 

2 

Hull  volume  inconsistent 
with  displacement 

PRELIMB 

3 

Longitudinal  center  of 

buoyancy  /  long,  center  of 
gravity 

PRELIMC 

4 

Error  in  range  or  increment 
of  variable  conditions 

PRELIMC 

5 

TDP  calculation  incomplete 

PRELIMC 

6 

TDP  file  lable  ^  title  data, 
col.  1-30 

Errors  in  the  calculation  of  the  two-dimensional  properties 
will  be  self  explanatory.  However,  if  an  error  is  found  in  the 
energy  balance  check  on  the  results  of  the  two-dimensional  lateral 
motion  calculation  the  message  is  printed,  but  computations  proceed. 
It  has  usually  been  found  that  such  errors  in  the  energy  balance 
check  have  little  influence  on  the  calculated  two-dimensional 
properties. 
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APPENDIX  A  PROGRAM  DESCRIPTION 

The  SCORES  program,  written  in  FORTRAN  IV  (RUN  Fortran 
Version  2  under  SCOPE  Version  3  for  CDC  6600  computer) ,  is 
structured  in  a  fairly  conventional  manner.  The  main  program 
serves  as  a  control  for  the  job  processing,  calling  various 
subroutines  as  required.  The  major  program  loops  over  ship  speed, 
wave  angle  and  wave  frequency  are  established  in  the  main  program. 
Data  are  transferred  among  subroutines  via  labeled  common  blocks, 
each  subroutine  accessing  those  blocks  specifically  required.  A 
special  common  block  labeled  PROGRAM  is  used  and  shared  by  many 
subroutines  for  storage  of  intermediate  calculation  data. 

Subroutine  PRELIMB  reads,  processes  and  stores  the  basic 
input  data.  Preliminary  calculations  are  performed  and  the  data 
are  checked  to  some  extent  for  self-consistency.  Subroutine  PRELIMC 
reads,  stores  and  processes  the  conditional  input  data.  Preliminary 
calculations  are  performed  including  spectral  density  calculations 
and  print  out  (via  Subroutine  PAR)  if  required.  Then  the  twor 
dimensional  properties  are  obtained,  either  read  from  file  or 
calculated  via  Subroutines  CKLEW,  ZIPSMO  and  TDLR. 

Subroutine  CKLEW  simply  calculates  the  two  Lewis  form 
parameters  for  each  section  and  checks  them  against  criteria  to 
insure  positive  contours.  If  necessary,  the  section  area  coef¬ 
ficient  is  increased  to  satisfy  the  criterion.  Subroutine  ZIPSMO 
calculates  the  two-dimensional  properties  for  vertical  oscillation, 
while  Subroutine  TDLR  does  the  same  for  the  lateral  and  rolling 
modes.  The  latter  routine  follows  both  the  method  and  the  notation 
of  Tasai.  Subroutine  MATPAC  is  used  by  ZIPSMO  for  solution  of 
simultaneous  equations. 

If  lateral  plane  computations  are  required,  Subroutine  ROLD 
is  used  to  calculate  the  natural  roll  frequency  and  the  additional 
roll  damping,  to  approximately  account  for  viscous  effects. 

The  basic  ship  response  calculations  at  a  given  condition 
are  performed  by  calling  Subroutines  ALINT,  COEFF,  EXCITE,  MOTION 
and  BENDSH,  sequentially.  Subroutine  ALINT  finds  and  stores  the 
value  of  each  required  two-dimensional  property  by  continued 
fraction  interpolation  in  frequency  parameter  (equal  to  circular 
frequency  of  encounter  squared  times  draft  divided  by  acceleration 
of  gravity) .  Subroutines  COEFF  and  EXCITE  calculate  the  coef¬ 
ficients  and  excitation,  respectively,  in  the  equations  of  motion, 
which  are  then  solved  in  Subroutine  MOTION.  Subroutine  BENDSH 
then  calculates  the  local  loadings  and  integrated  moments.  Closure 
results  are  calculated,  if  required.  Throughout  all  the  calcula¬ 
tions,  subprogram  function  SINT  is  used  as  a  simple  integrator. 

The  ship  responses  d  each  speed  and  wave  angle  are  printed 
out  by  Subroutine  TNIRPA,  including  closure  results  if  required. 

If  irregular  seas  are  used.  Subroutine  STATI  then  cc? ~ulates  and 
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prints  the  response  spectra  and  statistics  for  long  crested,  or 
uni-directional,  seas  at  the  particular  ship  speed  and  wave 
angle.  Only  the  integrated  spectral  response  at  each  wave  angle 
is  saved,  so  that  the  response  spectra  for  short  crested  seas 
are  not  available.  For  short  crested  seas  results,  Subroutine 
SPREAD  is  used  after  the  full  range  of  wave  angles  has  been 
depleted.  The  integrated  responses  over  wave  angle  are  computed 
and  printed. 

After  completion  of  the  specified  calculations,  control 
reverts  to  Subroutine  PRELIMC  for  additional  cases  with  the 
same  basic  data,  that  is,  the  same  ship.  If  no  additional 
computations  are  required,  normal  program  termination  occurs  in 
*  Subroutine  PRELIMC  upon  input  of  a  run  control  tag  less  than  0.0. 

Only  one  special  system  subroutine  is  included  in  the  program. 
This  is  referenced  in  the  main  program  by  CALL  DATE  (DTA ,  DTB) 
which  provides  the  current  date  in  the  argument  variables  as 
Hollerith  data  (DTA  =  MMMbDD ,b!9 ,DTB=YY) . 


Program  SCORES  -  Input  Data  Card  Summary 


P 

*3 

O 

O 

•H 

« 

K 

n 


<G 

p 

a 


TJ 

c 

o 

u 


Conditions 

Card  (see  legend  Format 


Number  below) 

Parameter.  Entered 

It  Columns 

I  1 

* 

Title  information 

A80 

2 

* 

Option  control  tags;  number  of 
segments 

1112 

3 

• 

Lengthj  density;  gravity;  displace¬ 
ment 

10X,  4F10 

4 

* 

Breadth;  area  coeff.;  draft; 
centroid  (each  station) 

4F10 

5 

OT  (E) >0 

VCG;  roll  gyradius  (ship) 

2F10 

6 

OT(B)-0 

Long,  gyradius;  LCG 

2F10 

7 

OT  (B)  >0 

Weight;  VCG;  roll  gyradius  (each 
atation) . 

3F10 

8 

* 

First  sta.;  last  sta.;  increment 
for  moment  calcs. 

3110 

9 

* 

Run  control  tag;  initial,  final  and 
incremnt  in  wavelength,  or  freq- 
quency;  initial,  final  and  incre¬ 
ment  in  speed 

7F10 

10 

OT  *:i)  >0 

Fraction  of  critical  roll  damping 

no 

11 

* 

Initial,  final  and  increment  in 
wave  angle 

3F10 

12 

OT(D)>0 

No.  of  spectra;  parameters.... 

110,  10F5 

OT(D)-3 

Additional  corresponding  parameters 

10X,  10F5 

Legend  for  conditions:  *  »  always  necessary  in  data  deck.. 

OT(-)<-  necessary  only  if  Option  Tag  indicated 
meets  condition  shown. 
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APPENDIX  B  -  FLOWCHARTS 


Flowcharts  follow  for  the  main 
program  and  each  subroutine .  The 
references  given  on  the  flowcharts, 
such  as  C-01  etc.  (above  and  on  the 
right  of  the  symbolic  outlines) 
correspond  to  numbered  comment  cards 
included  with  the  FORTRAN  source - 
program,  and  listed  in  the  next 
appendix. 


41 


rt,nwcw»0T  fop  M»jjN_'ponf,6iM_ 


>>>>>>>>>>>>>>>>>>♦ 

♦  ♦  "  *  C-o 


>>>>>>>>>>>>>>>>*<<<<<<<<<<<<<< 
-f --  -  -  i  T-flA  * 


”  ♦ 

♦ 

• 

*  SFT  ' 

4 

PFAO*  PPOCFSS 

AMO  STOOF 

♦ 

♦ 

• 

TOP  APPAY 

♦ 

RASIC  INPUT 

OATA 

♦ 

♦ 

Il14«r  t  twits  . 

4 

4 

♦ 

♦ 

4 

CALC  FPFO 

, 

ofpfnofnt 

PAPAMFTFPS 

*  c-oi 

•  • 

•  CALL  popvtmA* 


.•  CALC  •.  MO  ♦ 

MOMT  ALONG  ♦•>>>>>*  ♦ 

— ^71  - vt 


PFPFOOM  CAUCML AT  TONS  AT 

fatm  wavf  fpfoiifncv  ano 

STOPF  PFSIJITS 


PFAO  ANO  stoof 
rONOTTTOMAt  TMPIJT  OATA 


_ _ W<<4<<<<<<<<<<<  » 

-  -  ♦  ♦" 

f ••••••••••*•  4  ♦ 

•  • 

•  CAt  l  COFI  TMC  •  ♦  ♦ 

•  • 

••»••••«•••••  4  ♦ 

*  ^ 

_ 4  •  ♦ 

♦  V 

*••••••••••••  ♦  ♦ 

”•  FOP  LATFPAi  ■“  4  4 

_  »  P|  anf  calcs*  • _ •  ♦ 

♦  CAtl  onto  •  ♦  t 

•••••*•••••••  ♦  t 

4  “  ”>  ♦  ’ 

4 

♦  * 

_ _ *  r-o? _ ♦  ♦ 

•  ♦  * 

•  SFT  INITIAL^  . 

,  SMTP  SPFFO  .  ♦ 

. . . .  *  ! 

A  ♦  ♦ 

>>>>>>>>>>>>>4  ♦  ♦ 

♦  c-m  ♦  * 

••••*••••*•••**•«  ♦  ♦ 


•  4  4 

4YFS  4  ♦ 

4  4  ♦ 

_ /////////////  _  4  ♦ 

7*fl>  jNfoFffH  /  *'  ♦ 

/  OUTPUT  /  A  ♦ 

/  mfaotnc.  /  4  ♦ 

/////////////  *  ♦ 

4  4  ♦ 

4<<<<<<<<<<<<<<  ♦ 

T  *c-os  > 


StAPT  LOOP 
OVFR  WAVF 
FPFO  OANOF 


•  CALI  Al  TMT  •  • 

•  FOFFF,F*CT7F  * 

•  AMO  MOTION  • 


•  TOP  MO^FMT  • 

•  TALCS*  • 

*  TAI  I  PFMOSH  • 


SFT  INITIAL 

SMTP  SPFFO 


$s*jser_ 


4<<<<<<<<<<<<<< 

SO  ♦ 


J  *  •  •  x  •  •  9  *  •  t  •  J  •  • 
.  STAOT  LOOP 
.  OVFP  WAVF 
'♦  ANOI.F  OAMOF 


♦  ♦  ♦  YFS  ,•  M^PF  •. 

>>>>>♦>♦  <<<<<<•*  VA%,<  AMOJ.FS  ,• 


U0°F  •  .  YFS 

WAVF  FPFO  .*>>>>> 


OPTMT  PFSUI.TS  AT  Tm(S 
SOFFf)  AMO  WAVF  AMOLF 


•  • 

*  TAI.I  TMJPPA 


4<<<<<<<<<<<<<<< 
4  ♦ 

4  C-0Q  ♦ 


JwrPFMFMT 
SMT°  SPFFO 


YFS  ,♦  MOPF  •,  MO  *  ♦  *  • 

<«<<<•,  **,MTP  SPFFOS  ,•>>>>>♦  ♦<<<<<<<<<•  TAIL  SOO^AO  • 


MO  ,• 

<<<<<<<<<•. 


ShoPT- 

tpfstfo 

CPAS 


4  »'0 

♦<<<<<<<<<• 


TP- 

OFfilll  AO 
SF  AS 


**  fll  Ojltt  AMO  PPtMT 
^MOOT-TOrqTFO  CFA< 
PFStn  TS  at  Tm|s  <rrro 


TA1.CU1  ATF  ANO  OPTMT 
STATISTICAL  OFSUI  TS  at 
TUT*  SPFPO  AMO  WAVF 
AMO|  F 


♦<<<<<<<<<•  cam  STATf 


42 


44-4-4-4444444  4  V 


IA  - 

u;  * 


4 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


•  Z 

•  < 


o 
Z 
4  ■ 


•  Z 

•x 

D 

h* 

•  Ul 

cr 


V  V  V 

V  s  \ 

V  •£  X  V 

VON  \ 

V  I  \  h  X 

V  O  X  2  OC  X 

v  s*-*oo:  \ 

4  +  Nttl  UN  - 

sat  oi\ 
N  UIIN 
X  3\ 
V  ZN 
\  N 

V  N 
N  \ 
N 


D  ■*  4  ♦  4  4 
O 


a 

o 


a 

Ui 


A 

4 • 

<r 

4 - 

4- 

4- 

4- 

4- 

4  4- 

4-4  4 

4  4 

4  4 

4 

4  4 

4 

4  4 

4 

4  «- 

4 

4  4 

-4 

4 

■4 

-4 

•4  -4 

-4-4-4 

-4  -4 

-4  V 

A 

A 

♦  V 

A 

A 

V 

A 

A 

V 

A 

C 

A 

V 

A 

Z 

A 

V 

A 

• 

•  •  • 

•  • 

•  V 

\ 

•  •  • 

• 

• 

• 

•  V 

X  N 

♦ 

• 

• 

•  V 

X 

V 

• 

r*> 

•  -I 

*-  <x 

•  V 

4- 

V 

X 

in 

•■  4-  U) 

• 

o 

•  U>  4 

z  <-< 

•  V 

O 

V 

O 

X 

© 

•  tr  uj 

• 

• 

* 

V>  1 

»  H  Z 

UJ  4- 

•  V 

1 

N  K  lu 

U  X 

i 

U>  Z  D 

a 

1- 

• 

UJ  o 

*  <X  -• 

7  a 

•  V 

CJ 

X  Z 

Ul 

►*  X 

o 

K-  »~«  _J 

• 

</)  M 

D 

>- 

•  -J  C 

O  Uj 

•  V 

x  ►-  i r. 

<x  x 

J  4  «I 

4 

• 

u>  a 

a. 

• 

-4  -4 

•  3D 

1  z 

•  -4 

•4 

x  a 

UJ 

D  X 

•4 

4  4 

©  > 

A 

• 

< 

h 

z 

•UK 

*— » 

• 

x  Q.  o 

X 

it  « 

A 

* 

m 

►*« 

• 

U“ 

• 

V 

c 

O  X 

u  »- 

A 

« 

M 

4 

<  <  e 

• 

V 

a 

*-*  X 

UJ  •  D 

A 

• 

c 

•  o  z 

<  O 

» 

X 

a 

U>  X 

I  -j  a 

A 

• 

•  D 

• 

X 

't  X 

o  o  z 

A 

♦ 

• 

• 

• 

X 

■t  X 

>  •-• 

A 

• 

• 

• 

• 

X  X 

A  ♦  ♦  4-  f  f  ♦  4 


♦  4  f  4  4  4 


4  4-  4 
7.  li 

uj  t- 
t  A 


4  4  4  4  4  4  4  4 


■1 

'1 

J} 

i 

} 

) 


j 


| 

I 

jfj 


IRfflS?®S8roS>iS3^ffiK®ta(S^«BSS55B¥^% 


45 


S\ 

X  X 
*-«  X  X 
ox  X 

I  X»-«  X 
O  X*  2  «UX 
X*»  **  tt  X 


x«r  z  < x  ♦  ♦  ■*  *us  « 


Ifl  A 
UJ  A 
>  A 


»  4  4  4  ■*  • 


a  h 
»-  4 
4  D 

z  a 


C  ►* 

a  x 
»u  a 
N  UJ 


♦  ♦  4  4  4  4V 


C  A 
Z  A 


2  2 

•  Oul  ■« 

•»  *  a  ^  a 


•  u.  •  o  a 

>  •  Wl  Ul 

•  4  4  4  •  UJ  M 


•  X  l 

•  a  * 

•  o  ► 

•  X  * 

•  a  . 


•  IUD  >  -J**  • 


•  a  u  jo  • 


ja 

4  c 

c  u 


O  c 

V  X 
CH 


•  a  © 

•  tf  «u  4 

•  t  _• 

•  UJ  *-  u. 

•  X  UJ 

•  Utf) 


•  Uh 

•o>  J 

•  o  » tn 

«C  Uj 

•  *«  o  a 


•2  4  C 

•  H  X 

•  4  ,/>  •■• 

•  Hina 

•  S4l 
•CI4 

•  c 


<4444444444V44>< 


♦  ♦♦  ♦ 


♦  ♦♦♦♦♦♦4  V4< 


44444  444  4 


V  A 
U‘  A 
>  A 


V  A 
W  A 
>  A 


•  J<  • 

•  4  7tr  • 

•  **  c  o  • 


V  C  »fl  • 


•  o  o  • 

VO  «C  UlUi  • 

v  •  ja  o  • 

►  4  4  •  U  2  • 


U  O 

a  u 
oa 

7  U 


4  4  4  •  <  ► 


•  U  C  * 

«  v  c  : 

•  »-  u 


•  «/<»• 


i  «■  i 

u  u 

0  V  I 


•  :  5 

•  II  b  O 

•  u  o  c 

•  ooi- 


•  IT  Z  • 

1*2  0  • 

e  •  o  • 

i  •  »-•*-  a  • 
o  •  t-  o  =>  • 

•  UilO  * 

•  ZlflH  •  - 

•  Z)  2  • 

•  li  OO  • 

•  z  u  • 

•  o  o  * 

•  JO  * 

•  <  ac  * 

•  o  *  * 


o*  .01/11  • 
o  •  I  Hlfl  * 

I  •  X  z  • 

o  •  o  u>  x  * 

•  <S)  »-*  o  * 

*  -*  •  .on  •  - 

•  a.  • 

•  u.  a  • 

•  u  u.  n  • 

•  j  iu  * 

•  <  o  »-  * 

•  UUUI  • 

•  1/  • 


o  • 

*•  *  a  o 

i  •  OUJ 

u  O'!  u) 

•  JU.  o 

*  •  z 

•  »-  T  <t 

•  rt  ui  <t 


-  •  s  * 

-H  *  <4  *-«  </l 

I  •  »J  4-  2 

o  •  x  z  o 

%  V-  til  *-• 


o  •  «o  • 

*U3ul  • 

»  *  ■  -I  z  **  ♦ 

•  4  <J  X  • 

•  SJ  Ui  • 

•  3  10  • 


•  to  • 

rt  •  z  • 

.-«  *  til  o  • 

•  *>_!*-•• 

o  *  t  <J  *-  • 

•  4  a  <  • 

4  *  •  ©  1  •  ♦ 

•  UuJ  H  . 

•  _>>--<  • 

•  **  Z  O  9 

•  <J  *-*  X  • 

•  a  • 


o  •  i/i  • 

i  •  u.  ►-  * 

U  •  U  3  2  • 

•  a  a  <  • 

•  *-  h  • 

•  h  a  to  • 

•  ui  a  z  • 

•  to  >  o  • 


A 

V'* 

V/A 

U7A 

\v\ 

\  O'  <ki\ 
\o  o\ 
\c/'  ^ 


♦  V  1 


44444+4*4*  ♦  + 


I  • 

U  •  to 

V  ►* 
•  U  —3 

*  «  ul  O 

a  •  r  /) 

A  O  u.’ 


a  *  x  'n  • 

e  o  ^  %  •  V 

a  i  sh  >*  •  « 

►  a  u  v  *  x  aj\  •  J 

Uj  x  ~  c  o\  •*-•*■ 

♦  *4xxt<in,4 

x  a  x  «/*  x  •  ^  £ 

.  x  u.  v>  x  •  ® 

X  aJX  •  X 

X  XV  •  * 


•  •  *t/»  • 

•  c  to*  *ft  •  h  >  r*  • 

•  Q-  >-  ?  •  R  *  J  <1  >"  * 

•  uic  c  *  *  *  ;  a  o  • 

•  rsi  «  •  O  *  W  jo  • 

'•  .U»-*  *U  4  7  • 

•  ♦  +  +  *  *  X  Al  •* 

•  <*  3  J  •  A  •  O.4.* 

•  «U!3  •  A  •  U’  O  • 

•oxo*  A  •*-£• 

•  UJ  U.  -J  •  a  *5,"  » 

•  10  U  •  A  •«/>•-*  • 


A  4  4  4  4  4  444+44444 


«.<tt-t  ♦*♦*♦■*♦■*♦ 


4  *  4  4  44  4  +  +  444*  +  ' 


fclv4r*->j 


;  iLhfir'/ufeKv 


52 


APPENDIX  C  -  LISTING 

The  complete  FORTRAN  IV  source  deck 
listing  for  Program  SCORES  is  given. 
The  numbered  comjnent  cards ,  such  as 
C-01  etc.,  are  cross-referenced  on 
the  flowcharts  in  the  previous  section 


MPftftPAM  SCOWES  HNRUT.OuTEi>T.TAPE»aINPUT.T4P£8aOUTPUT,taPC10) 


COMMON  /  /  IPRtJI.IS.lO) 

COMMON  /  C0N0A  /  RI.RAMMA.RRAv.RO 

Common  /  **DT  /  *0A<14) iOTa.DTR.IR.IC. 10. It. !'. la. I*iII,IJ.STS<Sl 
rOMxnh  /  8*50*  /  RRLtOISRLtTMASStYNCRTtRSTARtll) .AREA  121) < 

l  SECOClll>.ORAETt*l)»ZiAR|fl)  •  Xt<l|iil||OiJ|M 

i  0wettH(2l),0MASSt2l) .ZmT(21).8RL(21) •ZCO.XNfRT. 

I  XZ»eBr»8MtMtXKRl»MAKKR]«!NC*CB»*0LnM 

common  /  CAsrtA  /  NN.0*M<S1)  '»VLt*l>*OMNE<Sl)  ,VMIn.v*AI,0CLV. 

I  nmA,wA0(25)  .NAN4I .MANRA.DVAUa.NMl , wO 1 20l  t*Lllfll 

COMMON  /  10IH  /  MC»veNtA*S(»l»10>*Kl.»KU»!0iItf 

CO***QN  /  MlMO  /  } AtNItOXl « V*MAM8*0MC(A«tf AVCN*CW*OIl |Z1*9) »EAc*MA 
OAT*  STS/IHANi  tO,AMR.M.s..A*AVR.  .AhSIR.  .8HAvl/l 0  / 

c 

C  ••  ••  tMEClAL  STSTfM  SUBROUTINE  vMlC*  RETURNS  CURRENT  OATC  ••  •• 

CALI  DAtE  irTAtOTK) 

C«0I  REAO,  PROCESS  AnO  f *0RE  INPUT  OaTA 
CALL  RRELIM8 
Sn  CALI  MfllNr 

1*  «  1E.RT.0  )  CALL  ROLO 

C-0?  INITIALIZE  SHJR  SRIEO 
V  •  VM1N 

C-OS  I  OOP  OVER  *Av£  ANRtC  RANOE 
AO  ^0  t«  «H«A 

MAMA  ■  NlD(!il*»I/l«0.t 

C»OA  Aft  TOt  ARRAY  Ut*8E  LIMITS 
ML  •  1 

1R<  IE  .OT.  1  )  ML  *  3 
MU  •  10 

IR  <  IE  .LT.  1  .OR.  AMOO’MAD<!W)illOtO>»eO,0*0  )  MU  ■  2 
t E  l  IM.LT.E  .or,  MA»AR:.C0.NJNMR!  J  SO  To  to 
RRl*T  t20.MDA.OTA.OTa 
PRINT  '?!«  V.«AOII») 

IE  I  I 1 .EC. 1  >  »RJN*  S2A 
•R|NT  t?3 

e-OS  LOOP  OvrR  M»VE  RRE3UCNCY  RANSE 
TO  no  A*  TO* I  »NN 
OMrrt*  a  Ow«((0> 

RAVEN  •  0ME0A»0RE0A/*RAV 
MVL  f  T 0 1 •  f.OMRt/MAVEN 
RLLMO)  •  RVL  A 1 0 1  /IPL 

C-OA  CALCULATE  EREOUENCT  PARAMETERS 
CM  A  SRAV/OMEOA 
w?  •  «AVENP|C*-V*COS<MAN«n 
CMME(IO)  a  ME 
m/n  .  «C»vC/*RAv 


RROAR*m  scons  UNRUT. OUTPUT. TARElAlNRuT.TARtSoOUTRUT.TARtlOl 
ICONTJNUEDl 


C-OT  RfREftHM  CALCULATIONS  AT  EACm  ERCOUEnCY 
CALL  ALTNT 
CALL  COfEE 
CALL  EXCITE 
CALI  MOTION 

IE  |  1R.LT.?  )  80  TO  80 

CALL  RENOS* 

AO  CONTINUE 

C-OA  PRINT  OUT  RESULTS  EOR  TmJS  SREEO  ANO  WAVE  AnRLE 
CALL  TNJRRA 

IE  t  IO.EQ.O  I  00  TO  SO 
EAC  •  (ll.O/IOlSRL*RRLI»|.0tPn«UO|P«2 
CALl  STATI 
tO  CANT  fN'lf 

IE  |  JM.LT.I  )  80  TO  100 

CALI  SRREAO 

c-ot  INCREMENT  Shir  AREtO 
100  V  a  v*OELV 

IE  l  V.Lf.VMA*  .ANO.  VMJh.NE . vm*Z  )  AO  TO  *0 
«A  Tft  to 

t?0  EORMtf  (  1*1.  1 JAt.  A2.  31.  A| 0.  All 

t?l  E0»***T  (  t*03Rff0  ■  .  ER.A.  AS.  13m«AVE  ANRLt  •  • 

*  et.?*  ;i*  uft..  moment  results  ) 

0?)  pntM*r  I  1*0 j  2 1 1 .  SAHVEHT1CAL  REnO.MY.  L*TER*l  RENO.MT.  TORS 

IIONAl  MOMENT  /  22*  HAVE  EREO.  STATION  .  3  iJOMAmPL  I  .*J0t  RMASC 

*  >  /I 

SE*  ?AR**T  l  !*••  84*.  IThINON-OImENSIONALI  1 

ENO 


53 


lltMftUTI«C  FR£LtMC 

ro»M»fN  /  COM)*  /  Rl.CAWHA.CSAT.RO 

common  /  mhoT  /  *#OA  ( U)  .QTA.OTR.IS.  IC.  ID#  It .  IF#  ]«•  |N«  1 1 , 1  U.STS  (!) 
COMMON  ✓  CASOA  /  •FL#DIS7L.TNASS#YNeRT,lSTAR(2i).A»rA(tl># 

i  sccociti>*o»ArT<ti)«2tAffun«xx  tin  •moon  » 

i  OwtliH 111). OWAS5 (II ).2»f till *4»L Jill tZC4»xW€*r, 

I  l2RCRT.tN.#(INi(RI.MAXKR:»lNCC£S#ROLnFF 

COMMON  /  MlMO  /  lA<NS#OM#V.VAN««O*f:CAiVAVCM#C»«0XX|2l*S)<FAC'«A 
COMMON  /  FR04RAM  /  lT0»A4ec4J*),r(2l»,m<lt)*Wlfl! 

•I  •  3.1*14414 

t*  •  o 

c-oi  »FAn  c *nO  ffoccsS)  VAsic  In*ut  oata 

1  Rf AO  «01#*OA 

»r*n  «ni*  u. jc. iCtiD.lt. If #I«.Im#ii#ij#n 

*  m  w*I  A 

t 9  (  M.OT.IJ  I  AO  TO  «tl 

Ml  •  M 

f>0  2  U»*m 

I  STA|T»  •  J-0 .SO* I 1 .0* 1 A1 

•fif)  *03#  *FL#6AHNA,MAV,0ISP(_ 

•fAf>  *0*.  (RiTAR(l)  iSfCOCCD  #0«A^T  ( I)  »I3AK(I) 

T F  t  2ClR(2t,»T,0*0  .OR.  IC.LT.l  t  CO  TO  * 
on  t  ui.m 

a  •  <i.o»?.e*srcocan/c.o 

If  {  A  »4T .  O.CO  I  A  •  O.tO 
3  70*R(D  •  ORAFT { I)  *A 

*  T f  (  It.vT.l  >  00  TO  II 

rfao  *e*.  tcc.raosro 

II  If  < I A.OT ,0  )  AO  TO  10 

•r«n  *ntf  raootR.ck 

AO  TO  1] 

10  »Mrt  *A*f  (0»eiiH(|>  #ZVTU>*CRl<I>  #A#1«)«M| 

If  I  a*L(|).6T.0.0  .OK.  JF.LT.l  )  CO  TO  11 
00  *  I»1#N 

*  A*Ut>  •  KAOCMO 

11  Rr*n  *04 #  Ml**Kl,tUX«KI, INCUS 

f-0?  •kfltminakt  calculations  u*ON  CASIC  INKuT  oata 
40  .  CAmMA/GKAV 
0X1  •  4KL/N 

\f  (  1C.4T.0  >  00  TO  13 

TMASS  •  OISKL/OKAV 

■  til)  •  IKKL-(l-lA>*OKl»/I.O-C«L 
YNFKT  ■  TMASS*KAOCTKMAO«TK 
00  TO  IT 

13  on  10  I«|«» 

If  t  IC.CT.O  >  CO  TO  14 

owfioHdt  •  d«cxshid*orav 

14  r>MAS«(I»  •  0»£10Mll/{**AV«OXl> 

TF  (  JA.tO.O  >  CO  TO  IS 

tr  (  I.fO.l  I  OKASStl)  •  DMAS* (1 ) *1.0 
T f  t  I.fO.N  )  OKACSlM)  •  OmASS (Hl *1.0 
is  vm  »  nMAS«tn«ii-i) 

14  Mflt  •  DMAS*  1 1 1 *ZkT (II 

TMA4S  •  SINT ( lA.M.OUASSiOXX ) 

HISK|  •  TMASSMK4V 

KIM)  •  SI"T(!a,m,y,dx!!«Oxi/tnass 

rn  -  lKKL-»l-XA»*Oxl)/t.0-KXm 
IF  I  IC.LT.l  »  CO  TO  IT 
*79  •  SlMTMA(M,tf(OU) 

IT  INfRT  m  TMASS**A04KO*KA0*K0 
OO  1*  !•!#•* 

•Kf A f 1 )  ■  CSTAK(X)A0KAFT(X)«SCCOC;X) 

Till  •  A»fAll)*(t-l> 

Mil)  •  KSTAKCI)#»3/lt.0-AKlA(X)«ZiAK(I> 

*!tf)  •  xX  tl ) •OiT*(X*l) 

1C  tTSOrXi  •  xitt)*xX(X) 

COM  ■  STNT 1 1 A.m. AKCA*0X1  1  •RAMNA 

COL  ■•1INTItA»M.TtDin*0xI*CAM«A/C0lS«(iKL.(l*tA)*0xTl/I.C 
if  t  IC.CT.O  1  CM  •  SXNTCIa«m.»«OXI)*CAMMA/C01S*ZCK 

f7K»#T  «  0.0 

IF  I  IC.CO.O  I  CO  TO  20 

e-03  CALfi»L*TC  LONOITUOINAl  MASS  MOMfNT  Of  InCKTJA 

00  10  t«1«n 

1*  vi!)  ■  OMASSllloiISOlX) 

VNf*T  •  SXNT (IA.M.TiOXI) 

RAORyR  ■  SOKTITNCKT/TMASS) 

IF  (  IC.LT.l  1  CO  TO  20 
7mTC  ■  WZS/TMASS 
no  if  WI.n 
7mT  id  •  JMT 1 1 )  -7mTC 
It  »(tl  •  OMASS(l) 4ImT <Xl*XX (X) 

IJKF*T  .  SXNT ( lA.M.tf ,0l! ) 


C-OA  FKlNT  OuT  SASIC  OATA  IlNCLUOfNC  »CtULTS  Of  FKOCCSSlMCl 
10  RRJNt  *FO»MO*#OTAtOTC  . 

MINT  COl#  !A»!K«!C»!0«!C*IF#XC#<m.I!»XJ»M 

FKlNT  C*n 

FKIN*  *03#  CFL#CAMMA»OtSFL*CKAV 

FKlNT  CO*#  (STAID#  •STARITlfSeCOC<D#DK»FTl!>#ZCA»(!)#0»CXCH<D# 

X  7*1(11#  C*L (I) # 1 *1 #mj 

IF  {  IC.CO.C  I  FKlNT  C04#  CCL.KAOOvK 

IF  t  IC.CT.O  1  FAINT  COS#  zcc.radcko 

IF  t  IS, CO. 7  • anO »  MXNMRI.CO.MAXKKI  1  KK]mt  SOT#  MlNK*I 
FKlNT  433 

IF  (  IC.CT.O  1  MINT  COS#  mISFl 
KKtNT  CIO#  CCL #COfS 

IF  I  IC.CT.O  I  MINT  |04«  CCl»KADCtK 

IF  (  IC.CT.O  )  MINT  CSC#  CM 

f-OS  CMfrv  MTS.#  VOU)#*C#  L.C.S#  ACAINST  ojsflaccncnt#  l.c.s# 

IF  (  IC.ro. 0  >  90  TO  21 

IF  (  aCSImISFl-OISKli/OISKl  #CT.  0.02  1  CO  TO  MO 
21  IF  (  ACS'COrS-OlSKLl/DXSFL  ,CT.  0,02)  CO  TQ  *44 

msK^  ■  ci' is 

IF  ,  AB«(CX»CCt.)/SFL  .»T.  o.oos  )  CO  TO  UC 
KpTliKN 

C-04  FKKA4  ST0FS 
V4C  T(  «  (1*1 
S4«  Tl  •  IX#) 

vso  ti  ■  ;r*i 

•SI  mint  «40#  It 
STAR 

•  02  FORMAT  (  SfMOOFTlON  CONTROL  TASS  «ASi:OCfCMIJ/ 

I  FIX,  1013*  ISI#  1THN0.  OF  STATIONS  ■  *  IS  ) 

COS  FQKmat  (  SHPLCNOTH  ■  #Fl0.2.  Si#  ChOCNSITV  •  •  F 11. A/ 

I  C«  OISKL#  •  •  F 1 0 . 7 »  SI#  SNORAvITt  •  »  Fll.4  I 

C04  format  I  TSmOSTaTIOn  scan  aRCA  COCr,  ORAFT  Z->AR  MpX 
XSMT  21 TA  SYR.ROLL  /(  FT. 2#  AF10.4#  F12.4.  |F10.«  )) 

•on  format  (  shoos  •  •  ro.s#  si#  ishctraoiusikoll  ■  «  fs.3i 

•04  format  I  13H0LONC.  c.c.  •  #Fa.3#40M  (FMO.  OF  MXOSHIrS)  LONC#  0 
XYCAATUS  •  •  FT.1  » 

•  nr  FORMAT  (  AM#  FCHCAlCULATC  NOMfNTS  at  STATION  •  13  ) 

•  OC  FORMAT  I  l H# #  TAX,  AMOM  •  #  F4.3  1 

CCS  FORMAT  (  44* »  ISHOXSRL.lMTS.)  «  t  FI0.7  ) 

•in  format  (  13MOLONO,  C.R.  •  #  FC.3.40H  (fwD#  Of  HIOShIRS)  OISRL 

I* (VA| .)  •  «  F10,?  1 

•30  Format  I  ITMOCASIC  INRUT  OATA  ) 

•33  FORMAT  I  lANOOCRIvCO  RCSULTS  1 

SOI  FORMAT  |  13A4#  All 

•0?  FORm*T  I  1112  1 

•03  FORMAT  t  10X.  3F10.S#  F10. 3) 

•04  FORMAT  I  AFlO.At 
•04  FORMAT  |  31101 

•  20  FORMAT  t  1«4),  1SA4#  *2.  3X#  AlO.  All 

V*0  FORMAT  f  3MOSTOF  JN  SUSROUTINC  PfttLlMl.  CRROR  NO.  •  lit 
FNO 


SUCPOUTINC  FRCLIMC 

COMMON  f  CONOA  /  RDOAMMA.CRAV.RO 

Common  /  mhOT  /  HO A (14) #OTA vOTR# IS#  JC# 10* 1C*  tF« IC# IM« 1 1 < 1 J#ST3(9) 

COMMON  /  SASOA  /  SFL#01SFL#TMA9S#TNCRT*SSTAR(21)#ARFA(2D* 

i  SCCOC(ll) .DRAFTI2I) #Z*AR(2l) «Xl (2j) (XlSOitl ) • 

X  ORCXCMIZD .ONASS(Sl) #ZwT(21) #SRL(21 >  *ZCS#XNrRT# 

«  XZF£RT«OM,NlN#(Rl,MAXSRI«INCRCS.ROLnFF 

TOMmON  /  CASOA  /  NN»0MIMSl}»«VLlSl)#0KFei91)#vNXN(VMAx*0CLV( 

X  NFA#VAD(2S)«WANSI#MAN0A«OVANC#NMX,wO(2Ot#RLL(Sl) 

COMMON  /  HlMO  /  XA#NS#OX| •V#WANO#OMCCA#WAVCN#CW#Of X til #9) «FAC*VA 
TOMMON  /  TOR  /  M  #SlM|2l|,SSSSm)#Nr#0NT(29) 
common  /  /  TOR (21 #25# 10) 

COMMON  /  RROORAm  /  STORASC (*SA) ,RSA (10) . HOC (S) 

OATA  NF  /  IS  / 

DATA  OHT  f  0.0*  0.01#  0.03#  0.04#  0.10#  0»1S#  0#2l*  0.2S#  0.3*# 

X  0.49*  O.SS#  0.4T#  0,87.  1.01#  1.29#  1.59#  1.95#  2.*S# 

X  3.09,  3.C.  4.7.  S.c.  7.1#  S,T#  10.7  / 

C-Ol  RFAO  ANO  RRfNT  CONDITIONAL  INCUT  OATA  CAROS 
20  RfAP  •07#  4A»SWLtSVL#DCLVL*VNlN«VMAX#0CLV 
IF  (  VA-C.O  >  40,1,7? 

IT  PRINT  •?OtHf>A»OTA*OTS 
RRINT  93? 

RRINT  RCT.  WA#SRL*S»L#0£LRL.VN!N.VNAX.O£LV 
TF  f  IC.LT.l  I  10  TO  2C 
Rf AO  ROT.  ROLO*F 
PRINT  ROT#  ROLOmF 
?4  RCAO  ROT#  VANCt#«ANCA«0MANO 
FKlNT  ROT,  NANSf #WANSA#0mAN9 
IF  »  in.LT.I  )  90  TO  25 

•C*n  ROC#  N«I # (*D  I  ! ) #  Xal • 19) 


mint  tn*.  n«j • (»o{ I j  « l*i *iei 
I r  (  lO.Nt.3  i  to  '.9  *5 
•tin  m,  i  ilm x i  «:*n*2o> 

mint  9e9.  <«oti > «:*u#io» 

C-03  INEyT  OITA  EMO#  CHtCt 
»  II  ■  t 

ir  *  NmKKt.NC.HAiMi  .a no*  meres. Lt.o  i  t*  •  s 

T r  (  S«L.NC.I«t  • ANO.  3ELvL.Lt. 0.0  3  !*  •  3 

ir  (  v**m.Nf ,yMA«  .AND.  OtLV.Lt.0.0  >  !*  •  3 

tr  (  V*Ntt.NC.VlN«A  *ANO.  0v  AnO.  Lt .  0 . 9  )  t*  -  3 

\9  (  0  >  to  TO  *so 

c-os  initialize  iano  cmc*)  Internal  kasanctcks 

N  •  N$ 

nf t  .  i.o 

not  ■  t.o 

\9  •  NkI.OT.10  >  00  TO  *91 

NN  •  (Am-St-U /nCLNL'1.001 

tr  i  NN.0T.5I  3  00  TO  9S1 

1 1  I  IQ.iT.t  |  00  TO  30 

e-ot  •»cltm!na«v  calculations  ro«  imcoulak  rives 
no  *>  !«1*nh 

22  ■  S>L*CCL*L*(l"l) 

CALL  *»* 

JE  I  tr.LT.l  I  00  TO  32 
<  a  OO.OOl/fiaANO 
19  (  K.LT .2  )  K  •  2 
tr  (  K.0T.12  )  K  a  12 
OvANA  •  AO. 0/A 

-AN/ii  •  uo.o 

tr  (  mAnOA-wAnOI  .to.  0.0  )  00  TO  23 

HANAl  •  0,0 
AO  Tn  32 

23  manat  •  90. C 
AO  Tn  32 

C-0 7  ••EtlNlNAKr  CALCULATIONS  ro»  ACOUlAK  NAyfS 
3A  no  31  1*1 »NN 

31  n-H m  a  S0»Tl2.0*KI*OAAv/<5rL*0CLNL*lt-l»H 

32  —A  •  ImINOA— ANftJ)/0*ANO«1.00l 

jr  i  NVA.0T.2*  |  00  TO  OM 

Oft  33  IhI.NkA 

hAO  ( T  3  •  KANOmOHANOrit-i: 

33  roNTtNur 

C-00  CALCuLA" *  ThO-OInCNSIONal  SECTION  MOECPTttS 
C  AND  CONVENT  TO  OI*»CNSIONAL  RESULTS 
40  tr  «  IO.OT.O  I  «0  TO  SO 
PALI  C*Lt* 

tr  (  IE. or.)  )  30  TO  *2 

C-10  vMTTCAt  OSCILLATIONS 
PALL  ZI»S**0  tOCTl 
ri  a  •I4B0/A.0 
nn  t«  jai.N 
r AC  •  ri*BSTANlj|**2 
00  *5  t a 1 , NE 

43  rnrijii.ti  a  TO»(j*I*l>*rAC 
tr  |  IC.lT.l  »  00  TO  A3 


C*ll  LATralL  ANO  KULLINO  OSCILLATIONS 

4?  cm  tol»  tot»T> 

r*0  44  J*l»N 

OMftA  a  50*T (MSTir (J) / (2.0*0* AVI 1 
ir  f  »Aftft  .( c.  0.0  >  00  TO  At 

Mil*!  a  rO*ANEA(JI 
rs* |4)  a  RSA(J)/M09 
MAlAI  a  ftSA(J)»*STAr(J) 

•  3*14)  a  tjl  (3*  /MOO 
MAIT)  a  Mi  13)  #ASTA»(  J) 

*$!(•»  •  »SACT>/A|00 

•3 A i 4 )  a  031(3) 

•  3lf)0)  •  »5MM 
nn  44  la] »Nf 

on  44  ■■3«10 

44  Tn»(.|,I.K)  a  TOO  ( J,I.K)MSA|K) 

44  rONTTNUf 


C'I2  H*ITE  Tft*  AONAY  ON  rtLt  (TAretfl) 

43  art  Tr  «JA)  ihja  (!)  *1*1  *3) 

-•Itr  <]*|  (((TftMJ.J.K) «J«1.M) .tal.Nri.Kal.lO) 

;r  r  orr.to.o.o  .0*.  oor.to.o.n  >  it  •  a 


C-13  rrtNT 
4T  »»)NT 
Ml-t 

nn  «* 
3  T  *  • 


OUT  TrO-OINENSICNAL  SECTION  MOrtOTICS 
•2P.HOA.OTA.0T0 

9*7 

J»)«M 

J“0  *3f • 1 1 • t  A  I 


••INT  III,  ATI 

4t  ••tNT  999.  OHT(l).  (TOrtJ,l.Kl«K*2*10l* 

«  (  OHY ( t ) *  (TCr(J.t«K)tK*l*10)«Ja2*Nr) 

•ft  tn  3; 

C-OA  A(Aft  Tor  AMAf  mok  Met  turtle) 
so  ir  i  X0.0T.2  )  00  TO  91 

■rift  do)  iNocm •i«i*ii 

CO  «>  la) »9 

tr  (  MOCIt)  .Nt.  MOAII)  )  to  TO  SAt 
32  PONTtNjr 

•f  Aft  (|0)  I  ( tT!)A  l  J«  I  »Kl  tjal  »H)  .I"l »NE )  <K*1.!0) 
tr  (  10.C0.I  )  00  TO  AT 

31  M  a  3 

tr  I  tV.Nt.O  )  90  TO  930 

•mom 

0»02  no  "ir*  ro»  NCa  basic  inaut  data 

l  CALL  MtLlM 
no  Tn  2ft 

c-04  r t»n»  STors 
940  1*  a  4 
iSft  U  •  1>*1 
931  MINT  940*  t* 

tr  (  tl.tO.fr  }  MINT  9A 1 ,  HOC 

to  srnr 

•32  rnruiT  <  /smocondItional  mruT  data  card  mint  out  /) 

90?  rO»MlT  (  9*10.0 
90S  roPMiT  (  no.  lors.l  ) 

909  ro»HiT  (  10*.  1 8*5. 1  I 

920  ro—iT  (  1  h  1 .  1 3 At.  A2«  3*.  MO.  A2) 

940  rftAMiT  I  39N0STftr  IN  SU3K0UTINC  MtLlHC.  CMOR  NO.  •  131 

941  rn*«AT  (  i9Horor  riLt  laicl  •  sx«  sit) 

997  fftOMiT  |  IMO.iOiOAHTaO-OlNtNSlONAL  SECTION  rKOACKrtCS  /4*»3mME0. 
*/4X,l27HrABAH,  A-MIMCI33I  A  (tAR)SQ.  M-VJ«  (S‘  N-Su4lSt  « 
*(3.4-1)  N||,»Hl)  t-SUAM)  N-SUtlB)  r-sUAC*.Sl  N-SUt(r. 

*91  ) 

99A  rftO-iT  (  4H  ST*  •  rs.i  ) 

999  rOAMit  |  flft.4*  1 JH  INf 1NJTV  .  9212.4/  I  *10.4.  10C12.A)  ) 

f  NO 


AytrnjTtNt  rift 

CftMMAN  /  CONOA  /  rt.0AMHA.0AAV.ro 

ftOHMftN  /  MHOt  /  M0AI14)  .OTA,DTN,lr.!C.IO.IE.IE.tft.tH.II,U.ST»lS) 
CO  — ON  /  CASOA  /  NN.ONm t31).*VLt51)»0H*ElSl).VHlN.VHAli0ELY* 

I  N*A.MA0(2S) .HANOI »9AN0A.0aAN0»N*t »¥v if 01 »«LL<*1 J 

CO— ON  /  3TAT  /  SreCN(10f31>«MO<0*10.2S> 

COHMftN  /  MOOAAH  /  STOrAOrt3«0)tVl9i)«KVST(10.S) 

OlMNStON  srctl?) 

DATA  SAc/tHNCuNAN.tHN  (195. **3)  (HA. OKIE)  . 

X  *MrjlrSO»»HN-MOSK.tHOM!TZ  •tH(lttA). 

3  AMTvO  r A .tHA AHtTC . tHA .  JSS.tnC  19tT/ 

nnouir  •  o»ay*o»av 

CONAT  •  O.OOOt2T«OSOUAr*rtAA3 

C-01  CALCULATt  HAVE  JI[CT*AL  OCNSJTy  *T  CaCh  ErtOUtNCT 
DO  AO  KKal.NN 
VfttTM  a  OHM (KK)aOMM (KK) 

OHSO  ■  OM(KR)AyOJTHAVO!TH 

c-02  ioor  over  HIND  sreto  ior  st*  sr*Tti  mnoC 

00  A*  tal.Nvt 

u  •  vDtpAl.tMOMtt 

no  Tft  |  10»  20*  30  >  .  10 

C-3-4  NfU-ANN  srtCTRUM  (1933)  (hAlE.  SO  ThAT  Sjf.  a  2  T J-M  K.h.S.) 

10  •errr  •  (-2.0aO»OuA»I/(VOITmauau» 

ArrcN(I.«P(l  •  (CONSTAt*rir09t»l)/(OHJ0a0HN(KTI) 
no  Tft  *9 

r-3-4  MCOAON-NOSKOEfTZ  SACCT*uh  (199AI  E0»  Elk.LT  MISEN  SEAS 
20  »OHEr  a  •« ?4* lOSOUAO/ (U*U*VO!Tm) ) a*2 

ArcCNd.KO  ■  .c901A0S0UAr«Esr (rover I/OHSO 

no  Tft  49 

C-j-4  tan  rA*iHtirr  srtetruH,  oased  on  sIOnIEicant  aivt  mt]0»*t  A no  hian 
c  HAVE  riTBIOO.  StNtLAA  TO  t.S.S.C.  NOMJNAL  d9tTl 
30  AA  >  0.?50*wDd  )*V0  C) 
r  •  :o*t 

M  a  (O.MTe*2,OArt/»0(Kl)Aa4 
MHfa  a  .0t/(VOtTM4vOtTM) 

Arfrad.MK)  a  AiatMtir(90MEr)/OHSQ 


4«  rnNTtNuf 
SO  CONTINUE 


t « f"  t  l-».\,tY-  * 


■Tl*JiV  MA  trv^f -- 


B- 


55 


(■OlKI.Ktoll.X) 

it*:«i(N»n 


•ts  tWfMTf  "AVI  iMCtlt  TO  OtTAIN  HAVC  AMBtltuOC  STATISTICS 
nCl  •  Om«i J)-0NmI|> 

no  *n 

no  9t 

99  VIUtSKCKKtbl 

•  91ntu.nn,v.ocu 

wvfT<«.?)  •  so*T(«yST(K«in 
«vtT<«.3l  •  tov*T<K.2>*1.21*9 
toy ST ( rt . 4)  a  •»tri«.IHT,C 
«vSTtK.9l  a  KiTlmtl'MMl 
AO  CONTINUE 

-0*  Mint  OUT  «4VE  $OtCT«i  A  NO  AMKUTUOE  !TtTl|f|{l 
aatNT  tlO.MOA.OTA.OTt 
JSAalO** 
lllatlf.) 

aaiNt  ioi.  lS»C<n*t«lSt*tSt*l> 

I r  1 10. lT, 3)  FAINT  104.  ftoft I A) »Kal .NMl I 
If  1 JO,‘_T.3 )  AO  TO  »1 
FttNT  lot.  (ao « I *xal .4»ll 
t  a  10*N«I 
MINT  10ti 
9l  MINT  107. 

MINT  Ut 
no  9f  Itol.NN 

'•(NT  lot.  pa« ( T  > *  I t*tC* 'K.tl.Ktol.NtoXI 
9f  PftNTtNur 
A A  J  NT  100 
On  AK  Ral.t 

tt  •* I  NT  103.  STS<K|.t»VfTtl.«l»ial.N»t> 

BftUAN 

100  fna««T  (  i h  1 

10|  ro*a*T  I  IMA.^i'tlHVtvC  J"tcT»AL  OCnSITt.  .«AA.tM  S»fCTAA/> 

tor  roa*»T  e  iifi*.3> 

103  rft»«uT  r  AX.  a«.  tofIT.*: 

lot  roa-AT  {  I T> «  ? AHAO*  *1-0  SafEOS  (SNOTS)  6?  /  It*.  lO'lt,)) 
lot  fOtMAT  1  «*,  TmsIO.hT,  .ft. 3.  *'12.31 
lOt  fO**»AT  (  «x.  TMMN.fE*.  .ft. 3.  **12.3) 

107  fO*M»-  (  1*H0  IKCTM  NO.  .  It.  til 2  1 

lot  roaMit  t  IJN  „AyE  mo.  > 

tt A  ro»-*7  (  INI,  1 lit .  42.  3*.  *10.  4I> 

fNO 


SU»anUTtNE  *OL0 

po— on  ✓  cOno*  /  a:.t*NM4»ta4v.ae 

COMMON  /  *4*04  /  *FLcDltfL.Ta4tt.VNeaT.tST4ai21l.4*f4(21l. 

«  *CC0e<2l».n«4rTi2l),2t*0<fl».lH2i).xlf0(2H. 

*  0vCI0Ni2n.0M4iS(2n.2aTl2n.»«l.(2l).ZCfl.*NetT. 

I  lZfe»T.tN,M!Nvai.N4IK*!.lNC«CS.»0L0ff 

CAMitrN  /  TO  I  ■  /  «e*«ZN.4N$itl.lO).KL.4U.10.1W 
rflaap..  /  MlK>  /  I4.NSi0*l.y.W4Nt.0N£t4.M4VCN.Ca.0l* (tl.tl «f4c*v4 
rO“**«N  /  aaooa»M  /  STOOAtEuSTi  .V(2li  *V(21) 

f-01  INITIALISE  a*4*aCTZ*»  atouiaCO 

*4,  •  7 
ail  ■  t 
tof  •  0.0 
•At  a  0,0 

C-02  Mlfltlt*  N*Tga*l.  »0U  taityCNCV.  INC!  U0JN4  40010  InE»TU 
2  tof*  «  nr 

tof  ■  toaTiOISFLtoON/iiNraT*at!i) 

If  <4tS(MEX/>c.-l.0>  »tT.  O.ol  )  00  TC  • 

tof N  •  toftotot/Oto*V 

£41.1  AUNT  *' 

00  *  1-l.NS 
VlT)  •  ANSI I  .7) 

*  toll!  ■  4NS< I »•» 

aA!  •  SJNT ( ! a.ns.y.O*! I 
B*D  •  51  NT  <14. NS. 4.0*!) 

no  Tn  t 


C-03  r*LCi.W*7C  ADDITIONAL  *OLL  OAMntNC 
10  afttnaf  •  i.n^OLOff  *0ISPL*tM/«r  -•*0 
•at*.  ta.  «r.toto*s,aocOfr 

tof*»>BN 

to*  ffiMMit  »  ///MX,  ttHNATUHAL  *OLb  f»COUZNCT  •  «f|0.*/  4*.  JAmCALCU 

II A  Tf 0  totvf  04toFX**0  in  »oll  •  .  tl*.*/  ItH  ACOITJCnaL  VISCOUS  04a* 

I JN*  IN  *OLL  •  .  f  I*  * *  * 

r»n 


•UNCTION  SINT  (INTn. JfT.OUl 

iNTfotATf  Th£  fuNCl.ON  Till.  VMICH  It  74lut*TC0  fO*  J  *0I*T|  AT 
fOUl. DISTANT  JNTE*VAIS  Of  0*1 

If  INTO  to  o.  USE  SlNfct  SUMMATION  TIMES  oil 
I*  INTO  •  1.  USE  TftAffZOlOAL  »ULf 

niMfNStO*  Till 
Su“A  a  9,0 

no  in  i*i. j 

10  tU*A  •  SUNA.Vt!) 

If  (  INTO  .ra.  1  I  SUNA  a  «UMA-,Y(l)«VMM/2.0 
SINT  •  n*!*9UN4 
•ttuaN 
f  NO 


SU*fnU71NE  CXLfto 

COMMON  /  CONO A  /  ft.t4MMA.ta4V.K0 

COMMON  /  OASOA  /  HfL.D!SfL.7MASS.TNEKT»tSTAa(2!>.AafA(tl). 

I  SECOE (21 > .OtAf T (21) .2(40(211 .21 (21 >  «*1S0 (21 1 < 

t  OMEIOM(2!» .nMASS(21».ZvT(2l».t*L(21J.2ct»*NtKT. 

t  XZK£«T*tN.MtNKK!.MAX«Kt.lNeKtS.KOLflKr 

common  /  TO*  /  M  »saN(2h *ftti(21).Nf.0MT(2S) 

COMMON  /  /  TOf (21 *2t. 10> 

OAT*  E»j.  E»0  /  tHlN  .  2hOC  / 

ataa  •  aia3.o/3>.o 

C-01  CMEfr  SECTION  • aKAMETEKS  AtAINST  lC«!S  f OKM  CKITfKZON 
no  1  i«i»M 

taaaditotccncdi 
tf  in*Af7(I>.LE«o.0>  to  TO  11 

*aHn)aasTAaii)/<2.t«oaAf  T(tii 
J* (AAM( I I.Lf.V.O)  00  TO  11 
If  <**rt(fl»i.O)  r.2.* 

C-02  7f»n  SECTION 
11  t*M||)  to  0.© 

00  Tr*  1 


•0  to  t 

tO  TO  t 


r  tf  (Attain  «iie.  fia*A(2.©*s*H(im 
(Midi  •  a!tfMto(?.o-S*Mtn  I 
no  -o  7 

0  If  I'MM!)  .« t.  »U»*(l.0.l,O/MMUIII 

s"*»in  •  ai«»*ir.o-i.o/s»Mim 
t  aaiNt  2A.i.f *i.*ecoc(i) .saiain 

no  to  i 

9  if  «  stand) .cf.  fI*(StM(n.l.o/Stnin*IO.O)/32.o)  to  to  i 

**A*<!)  •  • I* (S*h ( I > «1 .©/StHi It  *1 0.0) /32. 0 

aaiNr  to.  l.ean.sccoEdi.staadt 

1  COMTlNUf 
afTi-AN 

2©  fO»M*T  l*Mt  SECTION. I 3* AJM  NOT  VAL10  lC»IS  fOtM.  SECTION  **£4  CO£f 
If.  »A2.l2MC»EASfO  f aOa.f 7.*. Jm  TO.f 7.A.21H  fO*  T,0.a«Of.  CALCS.) 
f  NO 


*i»A*auT]Mt  ZIKSmO  I0ET) 

common  /  Toa  /  NS.SRH(21).Satt(21).Nf .0MTI2S) 

COMMON  /  /  TpA <21 .23.10) 

common  /  aaooMAM  /  sTOaxoi;  nt?» ,  00T. 

■  coiin.sKin.cofdn.sif  di).siRi(A.si.siKut«.«). 

f  S«t(  1 1 ). SM  111 )  *301(1 1 ). SSA  nil  *saA  dll. SOAdl). 

t  CfA<*.9)»eOAt9).Cft(9).tldO. 111. fax'  Sl.fOK  91 

OfT  .  1,0 

i  *rT  .if  TaitoNUMfTKic  functions 

no  i  1*1.11 

T|af»l 

COI 1 i aCOS (*l*0.l9707t> 

Md:»*INtX!»0.1tT0Tt) 

COf(I)toC0Sl*lto0. *712341 
1  n]A(f  ImsINOI*©.  4712J*) 

no  *  *•!.* 

K  to  **T 

A*  ■  AM*©. 1*707* 

no  4  !•) ’9 

Aftol-l 

A  *1*1 (f. 1 1 toSI N (  4M  toll. ©toAI.t. Oil 

no  t  j»i .* 

AJ*J 


% 

■1 

1 


•i 


V 

'I 


i 

s 


3 


i 


Bl  i. -.-r  1  AilK&a. -» LMCyMSaOk i*t£sj Vi.- 


T  S|«JIK.J)*SIN(  **  •!•••* Jl 

F  PONTtNJf 

C-U  LOO*  OVC»  *U«#ei  of  STATION* 
00  JO0S  K1*J#N4 
IF  <  Kl  .CO.  i  I  «0  TO  05 


f-oi  chccf  to*  comtaht  scction  mmhitmi 

*K  •  «M 

if  i  sM»(*i>.Nf .saabikkii  «o  to  ts 

IF  (  SBHIKl).Nf.SIM(«KI  I  #0  TO  BS 
00  O'*  IF  *»«NF 
TOFe«rl*lF«»l  •  TOFJFKiJF.il 
BO  TfiF<irI«lF.tt  ■  TOFIKK.JF.*) 

*0  »*  loos 

C-14  CMtCK  Fo»  Zf"0  scction 

OS  IF  (  SMBJKII  .lF.0.0  .0*.  SfHIRt)  .U.i 


.LL.0,0  >  SO  TO  BA 


SAN*  )• 141S4# (SBNB (Kl) *4,0-3.14 lSS|*StH(Kl) / 

) NSAM«Kl).i.O)*(S«M|Kl  1*1,0)  I 

IF  if  *A)  lt.U.Jt 
tl  MFITF  (4*1201  K* 

110  FOFmaTUTnOINCOFFCCT  FAFAMCTCRS.  ZJFtHO  OUITI  F0»  STATION  .IS) 

ocr  *  0.0 
40  TO  loos 

ir  <*!•?. mi«*se«mvA> 

SA*(StN|Kl)-1.0)/(S«H<Kll»1.0>*SAZ/SAN 
SI«UJ/|M*1«0 
SAA*SA*SA*S.0«S4 
S  AAA*fA*IAAO.  0*SA*fB 

4A7**(1.0*SAI*(0<SS)SS*0.04447*SA*O.0mTMAA«0.0lS4T«SAAA) 
1«4.0*SB* 10.2-0. l«fl4*SA-O.0ST0A«|AA-9.0141B*SAAAI 
SAN*. ( 1,0* SA|*(4, #44*7*0, •2BBT*SA40,01§OT4SA A >-O.0#fB 
l*(0. *4144*0. 03TflA*BA#0.01BlB»IAA) 
SFS*.a.O*SAI*IO.OtBST*0.01S47*SA)>B.04S44{O.OST04«e.llll44»A) 
fF4*.ll,0«SA)*0.0110T-0.0*S4*0.«141B 

C-OS  tOOF  OVC*  FUtODCNCT  BAN4C 
BB  no  1004  IF.I.NF 

SFWA  •  0«T(1F>*SBN(K1I 

IF  I  SFOFA  .ST .  0.0)  BO  TO  S 

SA#  •  0.0 

SC  ■  B.o 

40  TO  109S 

1  SN.SF4F t  / ( 1 « 0  *SA*  SB | 

SF1  ■  SA2-B»74f4l*Stf* (1*0*SA> 

SF2  •  san.o.tbs*o*s«mb 

SSBa.S.141BSMIn(SF4FA» 

no  l*  l»l.li 

•1*1-1 

St  •  B»*Ml.O*SA)4COII)*iF*COFII)> 

77  •  S«FMl.Q-SA)*Slft)-tOMIFII)l 
FMT«fAF(*VVl 

riacnsiiKi 

SI«.SJN(«*> 

SSB(T)*S.14IS»4«I**CHT 
SFBIfl*S,141***C t*CMY 
SOBm*SBtU)-SSB04(1.0*Bt*B.l> 

BA*B6BT<lini*7t*VV) 

BB1  •  BA 

JF  «ABf(TT)  »BT,  f.lt-t)  BO  TO  IS 

BC*1.BT0B 

BO  TO  U 

IS  FC*ASJN|1I/BA) 

14  40*0 
•f*0 

COBCF*COS(BCl 

COBCaCOBCl 

S!BC**B!*MBC) 

SjBCstJBCt 

JF  (BA  .BT.  4.01  BO  TO  IT 
AIM  *  2. B 

10  F04O0*Ml*CO0C2 
•f*Bt-BBk*fIBCr 
BB  J  HlB  l  *B  A*  ( AO**- 1 .  B )  /  A*«B ) 

COBCIvCOtClKOBCSJBCtMJBC 

StOCF*S!BCl*COOC»S!flC4COBC! 

C0BC»*C0BC1 

lnNatN*. 1.0 

jfiibbi)-(0.ooi*ba)i  n.ii.it 

tl  B04l-BO-ALO*(l.TBl4«A))nMT 
•C*I-BC*BCI4CNT 
BA**0*CI*B(4tJ« 

4S**f *Cl*BO*SJx 

Ao  to  ueo 

IT  BB1*|.0/M1 
00  IF  **•!»* 


BO*FO-COBCl*B*l 
•e*FC*SlBCl4BBl 
BB1*BB1*AMM/Ba 
COFei  *C0BCl4C0BC-StBCIMJBC 
siBCF*SlBCl4COBC*SJBC*COBCt 
fOBCF*COBCl 
11  CONTINOC 

AA*B0*).)41S014S14CHT4SIX 
AA*nrO.)4iSSi4Sl4llMV»Cx 
1000  SSA|J)*-BB 
S*A ( JI*-BA 
IS  CONTINUC 
SSA0*SSA(1) 

SFFUO 

SFOlaO 

SO*-A.OS1S« 

no  is  i*2.ii 

AJ*I 

SOA(t)*SSAm-SSA0*(1.0-IAl-1.0)/10»0l 

so*-so 

SF***C (1.0*SA)*SJ (J)*S.O*SBaSIF(J) )4|C»1STOTB*SO) 
«FQi*SFOl*SF0lI:*SFM 
SF*1*SF»1*SFA(J|«SFm 
IS  CONTJNUf 

SFQUSFOl-0  »S*S*4  (11)  4SFH 
SFF1*SFB1-0.S4SFA  1 1 1 1 *SFM 
F*A  I] .1 ) *SSAO 

roA.ii  *sseo 

no  ft  K*1(A 

MXN.l 

SlBaO 

Sl«aA 

SjAaO 

Sf  Aai» 

00  FA  Jal.S 
tS«»*l*J 

S> A*S1A*SC»< ISU*)*SIKI IK, J) 

SIF«S1B*S0B  ( I  SOB)  •SIKI(K.J) 

!•  FONTINOF 
00  F4  jal.4 
IS*  ‘  »2*J*1 

SFA.  ;4«S0B( !S0*)*SfKj(K« J) 

SlA*i  A*  SO  A  1 1  SUA )  *S  IK  J(K 
F*  CONTI NOC 

FOA  (*«>  •(  .144AAT4S1B*0.1S.T33SaS1B 

IT  F*A(N«,l>*0.1»44474StA*0.13SS3lASlA 

faAn.2ta.Sa 

FFA i f* 1 l •-! .O«0,1)3!14Sm 
F*A I  7.2) *-f A-O.Ollf l*St 
FAAIA.ll a-SAA-0 • 0040A*Sv 
F*A(«.f)a«SAAA-0.001S3*Sw 
F*A i| .S|*-0»3)S3S34S« 

F»AiF.S|**0.3B14T*ST 
FFA (4*3)*  1,0-0, IJ44t*Sa 
f *A |4*3)* -f A-0,4?3SB*Sa 
r*At«,*  *-SA A-0, 00BAI*SN 
FFA I l .4) *«0*tO*S* 
r*A,-»,4|*o.lsiKA*sa 

r*A|S,4)a0.1T»A**f« 


FAA IS. A » a-SA-0.010A0*BW 
F»A  M  »Sl*-0. l*Jt**a 
FFA(>,*)a0.Cf«03»S» 
f*A(S,9|*0»04TSF*Stf 
r*A|A.S)aO.H41T*S« 

FFA|A.Sia-1.0-0.0TAlB4S« 

FFBll )•) .0*SA*SB 

F*A  <F) *0,43441* (0,333333* (1 ,C*f A) -1 (BO*SB) 
r*B  111  aA.  31  A3  l*t  0,04447*0. 0A44T*SA*l,tMTl*fB* 
f **<« >«A. A34A1A (0, 00442*0 .00SS1*SA*0, l l 1 1 1*SB) 
F»A<S)ba.31F31*(O,OO7«3*O.O0T«34SA*O,OB1B2ASB) 
00  1  AO  1*1.9 
00  101  jal.S 
10|  •! II. Jl  a  CFA ( I , J) 

•  1  IT, 4)  a  -CQAITI 
00  1 0?  J*7,|0 
107  •) *|.Jl  •  -0,0 
>00  •!  (I. Ill  •  IF*| | ) 
nn  (of  1*4.10 
«  I-S 

MI  il,l)  ■  EQAIJJI 
on  1 04  jal.S 

104  Alir.J)  a  0.0 
no  |*>4  j»4*  l  o 
JJ  •  J“S 

10A  *1 It.Jt  •  CFA(tJ,jj) 

105  Bill, 111  a  0.9 


c-o*  C4ll  -*Tm4C 
C*U  *4TP4C 

tf  (OOT.EO.O.O)  OCT  ■  0.0 
00  J4  1  •  1 .1 

ppx(1»  •  lla»U) 

31  roiit)  •  •m*r,u> 

i  mpri-co*  1 1  j  -spoi  .cp*  (t>  *sm  *cpx  «3>  •in 

l*f»i  *R)MP3.CR*(5)«SP* 
4C«4Rf/(fl.T45**(1.0»S4.Stl*(1.0.S4MtJ) 
sA»*i.uM«*s-*«<}*T<e9xti)*e»xii)«CQxai*eoxmi 

C-PT  4T0Rf  RESULTS  IN  TOP  4RR4V 

1003  TftP(«lttftl»  •  SC 

1004  T0»!V|t|f «t>  •  S4P«S4R 
1001  CO^TlNUf 

RETURN 

twO 


SUMnvTlNC  n*TPaC 


C0»*QH  /  MOtOAH  /  STO«*Ptll1’).  0CT.JP4CE <101 ) .4 1 10 »l 1 > .SP4C* 1 10 > 
Of T«l . 

on  4  jai.t 
c-*"4(A(j,jn 

JPlaj.l 

no  1  MM9 

0al44(4(!,J)) 

Tfie.C)  4.1.5 
4  neta-OCT 

OO  T  A-J.ll 

4(f,*)a4(J,R) 

7  A(j,*)«n 

c»o 

4  CfNTINVE 

TRIAAtUlJ.jm  ZK.tO.11 
II  no  4  jajRj.io 

C0N4T4*(t«Jl/*U«J) 
no  *  K*j»i»it 

4  4ll«R)»A(!.R)«C9NST4A<J.R> 

IP (44S 14  < 10* 1 07) )  to. to, It 
It  oo  \>  lattlA 
*  ■  H-I 
4Pla«*| 

4-0. 

00  19  JaRPl.lO 
13  4a4*t ;K.J)*A  < J. 1 1 } 

It  4(4*11)  a  (4(4«1*.  >-S  1/414, 4| 
t?  RfTtlRN 

C-OI  PRINT  44PNIN0  »FSSAOC 
to  aaltr  «A«30l 
OfT  mO 
oo  TO  t» 

30  ROtaiT  (  34HOZOO  DETERMINANT  !»  SUM.  NATPAC  I 
en 0 


*uM00TJNC  TOUA  (OCT I 


e  tmt*  sur routine  h*to«ni  tmc  calculation  or  the  rotential  tmecrv 
c  AD9C0  NAtt  4 NO  WAVE  OANtjNt  AROAinTJIS  Of  TvO-OINfNIlONAL  LCUll 
c  fft4K«  IN  LATJ44L  ANO  »OLL  NOTION  NOOII.  Th|  MfTHOO  £N»LOTCO  It 
C  TMAT  Of  fUKUZO  TAfAl .  ♦HTOfOOTNANJC  fO*Ct  ANO  MONlNT  PROOUCtO  IT 
c  4V4V1NI  ANO  ROLLlNO  OSCILLATION  Of  CTLlNOZtl  ON  TNf  PREE  tv'MClP* 

C  IN  RrPORTS  Of  •C4C44CH  INSTITUTE  POP  APPLIED  MCCnAniCS*  RYUSmU 

C  HNlVMStTV,  JAPAN.  VCLUNC  I*.  NU*«l»  31.  1U1, 

C  4CI  ALSO  REPORT  if  J.  M,  VUtTS.  RTME  MvopOOtNAmJc  COtf^lCllNfS  POP 

C  4N4f i NO,  HCAYIN*  4 NO  ROLLlNt  CYLIN0CRS  IN  A  ERIE  Su*fACC««  REPORT 

C  NO.  |«4  (IN  CNOLtSH)  Of  LAB04ATOAIUN  VOOR  SCHCCMtOuNKUNOt. 

r  TfCMNlSCHl  MOtriCHOOL  OlLf T«  THE  NCTMCtLANOS,  J'NUARV  mt. 

f  ffMt.At*  1*70  -  0CC4NKS.  INC.  -  A,  !.  RAPP 

C  -  M3JCCT  NO.  ISO  (SSC-SAC  MOJfCT  SR-l74> 

c  444IC  INPUT  ANO  OUTPUT  VARIABLES 
C  HO  •  HALf-MCAOTH  TO  ORAET  RATIO 
P  4It  a  WCT1  ON  AREA  COCfflCltNT 

C  TI4  a  N0N*01M{NS ZONAL  P*EOUENCY  PARAMETER  (OMCOA-StUARiO  OVER 
C  BRAV1TT .  TINCS  HALf-tWAOTN) 

P  fill  •  AOOCO  MA|f  A NO  OAHffNt  4CSM.TANT  ARtAv  IN  NON-OlMCNSIONAL 
P  fOtM  (  A*  IN  VUtTS.  ABOVE  ) 


P  M  •  NO.  Of  T(ftM4  IN  P  ANO  0  (POLYNOMIAL)  SCtlCt  ( Jf  T  a  1] 

r  N  a  NO.  of  POINTS  ON  CONTOUR  fft  lCA»T  f QUART  f|T  (StT  a  IS) 

p  M|  a  NO.  Of  INTERVALS  POP  P.tLt-0  AND  X.SUf-R  INTEGRATION  (NjaNtl 1 

COMMON  /  TOP  /  NS*StM(21) .SPttlll) ,Nf .ONT(|N) 

COMMON  /  /  TOP ( || ,<3 . 14 ) 

COMMON  /  PPOtMAM  /  STOPatCl4*).AIlt.lO)*YM*).Tl  (U).PIt).S(t«l«H 

i  cotffint).  cocfpxiti.  seen  t).  tccti  ti.  pcoiuh 

I  P40(lt),Z(  tl.zil  f).XSll4).YSlltl«P|t),0(t) 

DIMENSION  CAM  1 141 

DAT!  CAM  /amnCnA‘1.  4HVE  CON.  4NT0UP  • 

X  4HUL-IE.  4HMAVE0  .  4MMATAU* 

r  AmAI  *  A,  4h3  CAlC.  Am  CPPOP. 

X  Ah<aEAaTI*  4HVC  fpf*  4M0UENCT, 

X  AMiHfPtT.  AM  tAL.  .  4HEPP0O  / 

P|-3.1A||ftT  f 

OfT  a  1,0  « 

II  a  0  < 

C-01  LOOP  OVfP  NuatfP  Of  STATIONS 
DO  JA*  41al.Nl 
MO  a  StM(Nl) 

410  a  SPPBI41) 

If  (  41  .CO.  I  *  40  TO  4S 

C-Ot  CmCCp  foa  CONSTANT  SECTION  PARAMETERS 
44  a  41-1 

If  I  SIRS (41 ) .Nt.SttS (44) )  40  TO  tS 

If  (  StM|4)).NC.SPM(44)  )  tO  TO  tS 
00  PA  If  al.Nf 
00  • A  JaS.lO 

•A  TOP(Vl.ir.J)  a  TOP (44. If . J) 

on  TO  109 

C-03  CMEOf  f OP  ZCPO  4ECTI0N 

•S  If  (  Sl8.PT, 0.0  .ANO.  M0.8T.O.O  )  tO  TO  tt 

OO  PA  Ifal.Nf  o 

no  PA  JaS.lO 
PA  TOP(vl.If.J)  a  c.‘ 

no  Tn  ins 

C-OA  COMPUTE  PCOMETk-lC  PARAMETERS  A.SUt-1  ANO  A.SUt-S 
At  It  a  l.A.HO 
XP  a  XAaXA 
XC  a  1,0-MO 
XN  a  XCaXC 
PR  a  XN/IP 

CCalA.RSIO*M0/(P!axp) 

AAaCf *40*3. 
apaR.ado.cC) 

CC  a  CC-4.0*mu/*4  t 

A3a(.OI.SOPT(aP4P|-A.aAARCC))/(E.RAA) 

A1  a  -XC*(1.T.A3)/XA 
A13a1.«Al*A3 
AA14.A13AA13 
TAJ  a  3. A*AJ 

C-01  CMfC»  ?MC  PCSULT 

IP  (A*S(MO-A1J/<1.0-AI.A)))  ,«T,  lO.E-t)  tO  TO  |« 

Tf  (Atl(SIt-PiaMOa(l,0-AlaAl-TASaAJ|/(A#oaAAl3l)  ,Lf.  tO.C-4> 

«  00  TO  30 

p-OA  fPPOa  RETURNS 
It  l<  a  t 
40  T»  a  ;«♦! 

aaiNT  47.  CRM(3alx-t).ERM(3alX-l)iEPN(3RlX) .HO.IIt.iIt 

DfT  a  0.0 
TX  a  0 

no  ta  ion 

C-OT  4f T  HP  VARIOUS  CONSTANT  fACTORS 
19  w  »  IS 

Nl  a  N  •  1 
EAC  a  Nl 

PN  a  Pt/(t.O«EAC> 

CC  a  PN/3. 0 
MM  I  P 

m  a  a 

wa  a  M.) 

C0N4T1  a  -TA3»»I/4,0 

C-04  CAlC.iLATf  t UNCTIONS  of  TMCTA  AROUND  SECTION  CONTOUR 
nn  4>  lal.NI 
PAC  a  I 
44  a  PNaPAC 
C44afOS(SS) 

CT4aCf)4  |T.a4S) 

SS4a4lN(SS> 

4T4a4IN(3,«SS) 

■  ft  a  ( (1 .*4 1 ) •44S-AT*ST5I /A  13 

Vrt  a  (ll.-4;)«CSS*A3*CTS)/4]3  .I 


IP  uftStxOI .LT.JO.E-A)  *0  >  0.0 

ip  um<»0)  .u.io.e-p>  to  •  o.o 
ir  1*0. lt,  o.o  ,ab.  vo.lt.  o.oi  oo  ro  «o 
Hit)  •  *0 
V3(»l  ■  ’0 

coerrim*  <u-41)*iss*ta3psts 

C^f/rl  ( J  >■  (l,*A3)*Al*StNtX.*SS)-J*»A3*SIN(4,4SS> 

C-04  CAiCULAU  P.SuP-0  COEPPICIENTS  POP  SftAY  4*0  POLL 
If  i  J.fO.M  )  00  TO  32 

4(1,1)  ■  VO 

4(  ! • vP)  •  XO*XO*VO*fO-1.0 

3?  continue 

c-io  L«n*>  ovr»  Focaur*CT  p4noe 
no  i*4  ippi.np 

!>•  ■  OMT(If)»Mn 
tr  |*IR)  4S,7g,M 
3J  rftH4T2«Pt**!»/»,0 

C-tl  CALCULATE  STREAM  4*0  POTENTIAL  PUNCT10NS 
00  4ft 

ro  -  XS(!I 
TO  •  rsm 

71  a  XO*XO*VO*VA 

Xft  m  /IP**0 
4(4  •  SIN(XK) 

ex*  s$  cos<x*t 
fry  •  exo<-in»YO> 

0-12  C41CIU.4TC  0  4*0  3  SERIES  POP  M4VE  IMEOPAL  APPROXIMATIONS 

ip  (vo. nr.  o.oooooo'.)  oo  to  33 
*1  •  pi/2.0 

AO  T*  34 

33  »I  •  *Ti*2(xO,Tft) 

34  M  •  xlft*30RT(XSI 
(A  •  X 4 

<C  »  ** 

*M  ■  l.ft 

00  •  0.3T72I3444*  •  ALOO(XA) 

•«  •  Xl 
es*  •  C03IXD 
M3  •  SIN(XI) 

CT4  •  CS3 
M3  .  333 
3ft  00  •  00*X4»f T3 
•3  •  PS.IA*3TS 
VM  •  XN*1.0 

IP  •  XMlC/XN 
14  •  XB/1N 

IP  II*. LT.  lO.Cf-T  )  00  TO  3T 

XI  •  C33«CTS-333«tT3 

3TS  •  3S3«CTS*C3S«3TS 

CT3  •  XI 

no  Tft  3* 

C-U  P*vr  ISTfOPAL  APPROXIMATION* 

JT  *4  ■PKV*(Q0*CX3*(PI-*I)«SXKi 
IP  •PKY»(QQ*3XK*(PS-P!)*CXK) 

C-I4  CftBPING  t£PP3  POP  MI  4 NO  PmI 
IX  •  XX4IJB 
PKV  ,  [KY»»1 

vil»  »  txr*cxp 

VUD  •  rPY#SXX  ♦  X4  -VO/XX 
*CO«t)  •  •CKYMIX 
P30IJ )  ■  CKV4CIP  -  XI  »XO/xx 
40  »’0NTINOC 

G-13  CO“PuTC  INTE0P4CS  POP  N. SUB-0  4N0  X,SUB*P  EVALUATES 
■4  •  PC0INl)*C0P7Pt<N!) 
wp  •  piooat*cftpppi«Nn 
if  •  PCft(NM*COPPPI«Nl) 
x*  .  PS0(N1)*C0PPP2(NI) 

*•  •  -1,0 
(.ft  *«i  I  ■!  «N 
#»  •  -P* 

•0  •  3«0*PP 

>4  .  X4«P0«PC0(I)*C0CPP|  I) 

IP  •  Xl*PO*P*Q ( I > •COCPP 1  ( I > 

•c  ■  «c.Po*PCon>*coiPHm 

in  •  *N.po»Psoa»ACOppp*(ii 
Yll)  •  Ylt)-Y(N)) 

4*  Y)  HI  •  Y 1  ( I  )-Y|  (N»  I 

C-l*  ftrTPpNIht  4LL  COtPPICICNTS  0*  p  4N0  0  SE*!CS 
i>0  *  !•!•* 

II.  »  I 
33  ■  PNAP4C 
44  •  *1,0 


»•  ■  COS (S3) 

ps  ■  2,ft*pa«P0-i  .o 

vP  •  PS 

*0  •  ?.0*ft$*BS«1.0 
It  •  SIN(SS) 

IX  •  2.0*X I*Pn 
IP  A  XI 

»  ■  2»0*XX*P3 

oo  ■  o.o 

00  30  JAMM.M 
00  »  00  •  2.0 
31*  ■  41/(00*2.0) 
fix  •  T43/(00*4.0) 

PP  r  YPAPR-XPAXI 
PP  ■  *0»>S>1IUI 

4(1. J)  ■  BB  •  (XlP/413) * ( YP/00.P0*SXK-PP*CXK*44* ( 1 .O/QQ-SXK-rXrt ) i 

VP  ■  PO 
PO  >  PP 
IP  ■  IK 

IK  m  XftPPS* YP»XX 

IP  (  I.NE.N  )  oo  TO  30 

PP  ■  00*00 

PO  •  (00*?.O)*(OQ*2.0) 

PP  p  (00*4.0)«(00*4.0) 
rt  p  (0o*i.ft)p(oa»i,0) 

3rci(j)  p  XIB»A4*((1.0/(EP-1.0)  -Al/(CQ-1«9)  -TA3/(rP-l .0) !• 
x  (1*0-41)  «TA3»(-1.A/(EP-B.0)  *41/(20-0.0)  *TA3/(Eft-V*0) ) ) 

3FC?(J)  p  *4*(2(ftPAlP(l,o*4))/(ES-4,0)  «B, 0*43/ (C3-14.0) ) 

30  44  ■  -44 
33  CONTfNue 

C-l?  30LVY  SIMULTANEOUS  EOUATIONS  POP  P  4N0  0  PtRtES. 
c  PftPM  MYI  COtPP ICIf NT  NATPIX  BY  LEAST  S0U4PES  METHOD 

*2  Oft  T  IpI.m 
ftft  T  JpI.m 
3(I.J)pO. 

00  *  3*1 *N 

4  *  1 1 • J)P3(I, Jt*A(ft.ll*A(K.J) 

T  3(J.I)PS(I.J) 

r-1?  PO®**  B.M.S.  M  vECTOP)  by  least  SQUARES  *»£TkCO 

00  4  IpI,** 

7 (I )pO. 

71(I)pO. 
no  a  JpI.n 
7(i».2(n *a(j,i,py(j) 

4  71 ( I )  •?! (|)*4IJ,1)*Y1(J) 

e-IT  INVPPT  COEPPICtrNT  MATRIX.  IT  REPLACES  0PI0IN4L  MAfPlX 
ftO  10  I»1*m 

ftft  «  JPl.M  A 

ft  ft(J,MP)«0,  * 

3*1.  **•)•!. 
f'lvpt'i.n 

IP  (4BS(0IV>.lT.  10.f-16)  00  TO  32 

Oft  P  JpI.np 

*  3(I.j)pS(I.J)/0IV 
00  7  JPl.M 
IPII.EO.jj  CO  To  2 
PACpft ( J, I ) 

ftft  ft  ftp) .MP 

ft  3(J,iOAft(jtP).S|I,K)*P4C 
7  CONTINUE 
Oft  P  J«),“ 

*  31  >.T)«S(J*MP) 

)0  frtWTtNUP 

P-l*  CALC'iLATF  P.SW3-2M  ANO  0.SU3-2M  SERIES 
no  )1  Ip1,m 

pin  >o. 

OlllpO. 
no  11  jp!*“ 

•  (i)>»(n*3ii,j)*2(j) 
ii  0(np0(n*S(i.j)p?i(j) 

r-l«  CALCULATE  N. Sup-0  .  m,SUB-0  .  x.SU8-R  4N0  -,sub*p 

PP  p  0 ,  A 

•ft  P  O.rt 

•p  p  ft.ft 

•3  p  O.ft 

ftft  Pa  jp>.«* 

pp  •  PP*P(j)«Sfri (J) 

PO  P  P0«0( J>«SPf I (j) 

PP  •  PP.P(J)*SEC?( J> 

PA  P3  P  PS*ft( J)*src7(J) 

PN**  p(-IA«CC*C0*STlAP(2n/4)3-PP/4A13 
rvft  ■(-ippCC*C0*ST1*0(2))/A13  -P0/4A13 
I*  pHC*CC*CftNSTfA(41»P(»,  -43*P(3n*PP)/AAn 
vi  piIN*CC*C0NST>*(4)P0(2)-43p0(3» )**S)/4413 
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C-19  PO***f  Nf  TC*MS  f0a  FINAL  PESW.TS 

»»  -  «o<ii«om 

oo  •  r-o  api!)  -fno  *9(1) 
pq  •  fno  1  >  *fmo  •uiii 

•  •  ■  I*  **U)  •  Y*  *0(11 

M  •  YR  ••III  •  I*  «QII) 

jr  (XM,P0,X)  90  to  at 

'•i)  »•*» 

•  I  •  A*S<OC*2.0/CPt*PII  -1 .0 » 

! t  (H.IT.0.Q2SA  |  90  TO  *1 

jr  (  A9SI  GO-4, 934#I/AmAa1 (A9StFK0*P(l )  » |A9S|PN0#0(1  J )  I  .AT.  O.lO) 

*  I*  •  I 

»*  a  hO/<SIO*FP) 

H  ( 1 1  a  ««*»0 
a»3j  a  X*»FA 
XK  •  XKa«GR7M|A) 

»l»l  a  4 

•<*>  a  >«*»S 

c-zi  snitp*  coefficient*  for  surseouent  roll  or  skat  c*tc- 

«a  on  AA  1*1 «N 
r<  a  acj.mPj 
All. at)  a  A C 1 • 1 ) 

AA  All.n  a  XK 

tr  (  MH.fTO.  XI  90  TO  H 
MM  a  t 

no  to  77 

e-22  POLL  RESULTS 

AT  IT  a  AAs<P$*»,0/<PI*P!)  *1.01 
ir  a  mO/(<.,0*5M»M> 

•  |A|  a  ik*FP 
■X?)  a  IKAFQ 

t«  a  XKaSQPTUfA) 

»|M  a  PI*FI**K/A,0 
a  I* i  a  VRaQQ 
00  TO  la 

e-zs  7 fro  frfouenct  calculation* 

70  7 A  a  1.0-Al-AJ 
»C  a  »*U 

•Cl I  a  »I*lCU9-*n*«taTA3«AS)/(*aO*HO*SIMxei 

IA  a  Al#M.O-AXaU.O*A3>4«a|a.0*3.0»»U/*.OJ*A3*«0.90*U.0aAr/7,) 
TC  a  I.OAMO/CSIn«AA13aAl3) 

•111  a  -iC**#/3.0 

•  |A>  a  lC*(<AXa(t,0,A)naa}.9,oaA3aUx«Ala(Al42,0))/9.9l/INIaAl)) 
*|T)  •  A(3» 

00  Y>  T-2,9,2 
T7  aCl  a  0,0 

no  Tft  JP3 

C-Za  r»»oa  RftuR*S 

99  IF  t  n.ro.o  )  oo  ro  103 

H  II  a  II  •  ? 

Till  «  II  *  I 

PRINT  97*  faMi3*Ix-*1,E«s(34ix-i>.t*N<3MX>.hC.S!9»*!« 

IF  I  IX.LT.*  }  90  TO  *3 

■  ■Tat  99*  P*  U 
aatNT  91 «  FNO.FMO.XR.YR 
93  I F  I  IX.LT.S  I  OCT  a  0*0 
II  a  0 

e-Z*  9T0*r  RESULTS  IN  COMMON  4**11 
103  mm  a  t 

00  A*  Tat. I 
J  •  1*1 

A9  TpMflfi.triJ)  a  9(1) 
xo*  roNTiNur 
IP*  PCNTTNUe 
i*f  Than 

97  FORmaT  (32MPJTOM  IN  SUMOUTINC  TOL*  OUC  TO  •  3*9  /  R9N  PARAMETER* 
I-  M(|  a  ,  no, A*  3li  9N*U  a  ,  F10.9*  3X.  9HXJM  •  .  Fit. A) 

99  FORmaT  |  UN  P  J!»I»  -  ,  9ei3.4  /  11M  0  SERIES  -  t  9CI3.A  | 

91  FO»“aT  IlM  NO  a,  tit.*.  AH,  HO  a*  112,4.  AH,  X*  a*  fJt.A.  A-.  YM  a 

x  •  eit.ai 
fNO 


SUBROUTINE  ALINT 

C  JN'rRPOLATE  ALL  "CQUl*CO  Tao-OJMCN»»ONAL  PROPERTIES  AT  PARTICULAR 
r  FREQUENCY*  70N  ALL  SECTIONS,  use  CONTtNUtO  FRACTION  NCTmOO*  «!Tm 
r  Six  POINTS,  TNAfC  ON  EACH  SfDC  Of  RIVEN  POINT.  ADAPTED  FROM 

e  SUPRPUTlNff  ACfl  ANO  A T|M  OF  *9VSTEM/390  SCIENTIFIC  SUBROUTINE 

r  FACRASE*  Vl»S,  !!!•(  ISM  Fun,  no*  NZO-OIOS-3  (1999)', 


COMMON  /  M1SDA  /  8FL*0lSFL.TMASS»YNCAT»9STA9l*li,A*FAItJ|* 

1  sccoEtzn  .o»AFTcin  *ziar<zx ,  ,xi  tm  •xxsotzi  i  ♦ 

X  ChEI0H(21)  ,0MASS«1)  *ZhT(IX)  ,9RL(Z1)  *ZC9*XNE^T, 

1  x2fc*t,om,minr«i,mai«r:, inches, rolopf 

common  /  TO**  /  NS*SSH(21 > *SRB9 til )  »NF*OmT(Z*> 

COMMON  /  /  T0PC21. 25,10) 

COMMON  /  TOIR  /  MC>aCN,ANS(Z;*10)«XL«KU,IO,IN 
OtMfNSION  YAL(i) *A*0(9) 

P-01  ton*  OVER  NUM9EA  OF  STATIONS 
00  2P  RlaltNS 
km  a  KL 

I  a  aEN*0**AFT  IK )  ) 

P-02  PmFPk  F0»  2FR0  SECT , ON 

IF  I  0»AFT,n>  ,«T,  0,0  )  00  TO  1 

00  7  KaKM.KO 
7  AN* (if  1  «K|  a  0,0 

no  to  Zft 

P-03  PMEPK  IF  X  IS  IN  RANOE 

1  IF  t  X  »LC»  OMT(NF)  )  00  TO  5 

2  JJ  a  NF-3 
00  TO  4 

C-04  MINAar  SEARCH  IN  OMT  ARRAY  FOR  X 

*  J  a  iiF 

»  a  1 

M  •  (J»l)/7 

IF  (  X  ,0T,  OMT (K )  >  00  TO  A 

.1  «  a 

AO  TO  9 

•  »  a  a 

9  IF  t  IARSIJ-I)  .97.  J  >  00  TO  9 

JJ  a  I 

C-OS  OMTjjj)  XS  JUST  tELOF  X  |N  ARRAY  OMT 
IF  (  JJO  .6T.  nF  )  00  TO  Z 
IF  I  JJ  .OE.  3  >  00  TO  4 

IF  I  JJ  ,0T.  X  .OR.  KL  *0T.  1  I  00  TO  3 
P-OA  FAFOUENCT  FanamfTCR  NEXT  TO  OMT (XI*  ZERO  FREQUENCE 

KM  •  ? 

TF  (  X  .OT.  0.0  >  00  to  31 

INS(kX.X)  a  I.0F7S 
nn  Trt  3 

33  *M  «  1,a*BSTAN(*1)/(2. 0*0* AFT (K X ) ) 

ANSCKl.il  a  TOFCRltZ*nAt0.23-ALOOtX*XHn/(O.Z3-ALOA«OMTt:5*XM)) 
1  JJ  a  3 

C-OT  SFT  ,jP  LOOPS  ANO  STORE  VALUES  FOR  INTERPOLATION 
4  00  1 A  KaKM.KU 

**  a  JJ-3  »4/|K*JJ> 
l»0  11  la]** 

TF  a  KK* I 

VALITI  a  TOP (Rl • lK.KI 
IX  AAOtT)  a  OHTCIKI 
•1  a  1 » A 
P>  a  VALCH 
01  •  0,0 
07  a  1." 
inn  •  1 ,(7* 

C-OA  PCWTtHyfO  FoiCTlON  INTERPOLATION  LOOP 
00  \»  TaZt 9 
JF  a  1-1 
00  IS  J*1 « JE 
m  a  VAUD-VALCJ) 

IF  (  M.NC.0.0  )  00  TO  U 

VAL(T)  a  1,P7> 

On  Tn  is 

14  VALtM  a  CAPOC  U -ARQ<  J|  )/M 
13  PONTTNOF 
«m  •  IHN 

•  1  a  vAl ( 1 > *PZ* ( X-ARO II-l))aPl 
01  a  V*L  1 1  >*02*  |X-APO(  1-1  n«0! 

IF  (  03, NF .0,0  )  00  TO  IS 

XMN  •  I, ITS 
AO  Tn  17 
I*  *MN  •  PJ/Q3 

IT  IF  {  A9SI1.0-XHM/XH)  ,LT,  0,0Z  )  00  TO  ZZ 

PI  a  PE 
P7  a  B3 
0)  a  02 
07  •  03 
12  CONI  t NI)F 

C-09  S700F  PfSULTS 

ZP  ANS(Fl«K)  a  XAN 


-r«TK5W^)  pm*-* 


\*  continue 
|0  CONTt^Uf 

Of  turn 

fun 


RuMnuriNt  cotrr 


COMMON  /  CCNOA  /  PI.JAMMA.8RAV.RO 

COMMON  /  tASOA  /  *RL.0:SPL.TMAIS.VNCRT.1STAR(21) .AREA  til). 

I  sccocttii.ORArmii.zURrtn.xmn.xlsoitlH 

i  ovtiOHjtn.OMAssrtii.zvTitii.iRLitiNZcitxNERT. 

*  XZFCRT.aH.NlNVRX.MAXKRl.lNCRCS.ROLnPF 

COMMON  /  TOJR  /  VE.VEN.ANSI21.10l.rL.rU.I0.lv 

COMMON  /  MIND  /  lA,NS.0xX.V,WAn8.O*£iA.vAYEN.CV.01*<21»S' .EAC*VA 
COMMON  /  COMO  /  CY(l2>*CLt2TI.ZV**»*TV«NV.KV 
COMPLEX  2V.MV.YV.NV.rV 

COMMON  /  MiOORAM  /  $TOIU0EU4tl,rmi«rxriO>»CX$IM,O'<5»»OCX<»!* 
x  orxsiu.Yrtn 

rt  •  l a 

u  m  mS 
nr  a  OXt 
TV  •  *.0*V 

C-01  CALCULATE  REQUIRED  XNTCORALS  OVER  SHIP  LCNOtH 
no  in  KaKLtKU 
no  1  UI.n 
7  vm  ■  ansu.m 
rrv)  •  xintcit.m.y.ox) 

tr  (  (K«1I/?,CQ.4  I  00  10  10 

no  A  lal.M 
4  v(!>  a  vmaxim 
ritll*  SINT(tT.M.Y.OX) 
tr  (  K.8T.4  >  00  TO  lo 

no  6  lal.M 
A  Y(!)  a  Y(!I«X1II1 

rxs(*i-»»iNTar*N.Y.oxi 

|0  CONTfNUr 

*1  a  SINT  (XT i“«MTAP*OX)*R*MM» 

00  10  lal.M 
is  Yin  •  ASYAaaiaxim 

RX!  a  SINTrtT.N.v.OX) *OAMMA 

no  14  lal.M 
14  veil  a  Y(l)*Xim 

PXSIa  SINT ( It .M.Y.OX) *#AMNA 

c-02  increase  roll  dampino  ito  account  roR  viscous  tEElerll 
k»i  a  r r«)*R0Lorr 
IF  (  KL.8T.X  I  00  TO  10 
fAC  a  «0/SQBT|VEN**3/aPAV> 

E|2»  a  f  |2)*FAC 
rxtX)  a  ex<2**eac 
r«*(»  a  rxsitiaFAC 

c-01  CALCULATE  PrOUlPTO  OCRlVATlvCS  ANO  TH£IP  InTCORALS 
11  TUX  .  t. 0*0X1 
MM  a  N-l 

00  ?0  KaKL.«U*2 
MV  a  (K»l>/* 

*«1X(1.KV)  a  (ANXll.K)-ANI(2,K|)/0Xl 
01XIM.KK)  a  l AN% (M- 1 *X1 - ANl (M i V 1 J/OX* 
no  XX  Tat«MM 

Zt  ntV(T«KK)  a  (ANR(l-liKI-AN|(X«l*K))/TOX 
on  X*  tat«M 
*4  Ylll  a  mxii.kki 

nrevvi  a  iint (IT.m.y.oxj 

tr  r  xR.ta.A  >  oo  to  to 

on  *«  lal.M 

x.  vet»  a  vm*xirn 

nnrvKia  sintmt.m.y.dxi 
tr  I  KX.OT.t  I  NO  TO  tO 
nn  XA  !■!«** 
to  Yet)  a  V  I  X ) *X l 1 1 ) 

nriR |XK) a  sjntut.m.v»dx> 

*0  fnMTlNUC 

tr  I  KL.0T.2  I  90  To  30 

C-OA  ro*M  COfrricIlNTS  roR  VERTICAL  RLANE  NOTIONS  IHtAVE  *  RITCN) 
cvr  1)  a  TMASS*r*l) 

cvr  xi  a  rexi-voni) 

CV<  s)  a  SI 

cvr  *)  a  rxeii 

cvr  41  a  rx(t)-v*orx.i)-tv*r m 
rv r  a)  a  ixi-v*cv r t) 
cvr  Tl  a  YNFRT.EXS 1 1 1 

.  cvr  at  a  rxst*)»v*Drxsri)-Tv*rxtl) 
cvr  ai  a  Bxsi-v*rx(ti*v*v«orxtii 
cvfio)  a  run 


cvnn  a  rx(ti-v*orim 

cyuxi  •  **i 

If  <  VU.lT.3  I  00  TO  40 

C-05  rnRM  cocrrs.  roR  LATERAL  RLANE  MOTIONS  iSVAY.  tan  •  ROLL) 

30  C*.(  1)  *  TNASS*E  r 3) 

cir  x>  ■  rr4>«v*orrt) 

CL  I  11  a  0.7 
cl r  a)  a  rxrM 

cl i  R>  •  rxi4)-vRorxr*i-TV*r'3) 

CL  I  A)  a  -V*CHt) 

CLI  T)  a  -r rR)'7C0*r<3) 

CLI  *>  •  -Fll0l-ZC8*CLl2|*V*0ElJ> 

CL  I  «)  •  0*0 

c<  no)  •  rxi3> 

clui)  ■  rxiAi-vrorxil) 

CL'IRI  •  0*0 

CLUI)  a  YNEAT«rXSI3> 

cl < i a)  a  rx$i4).v*crxsit)*Tv*ri<3) 

Cl ( 1*1  a  -V*CL 1 1 1 1 

CLUO)  a  -XZRERT-rX'f|-ZC0*rx(3) 

CL  I IT)  a  -rxtl0)-ZCQ*Cl(ll)*V*0rXl9) 

CL<1R»  •  6*° 

CLUO)  a  -rtSI-7CO*rt3l 
CLIXO)  a  »F IO)«7C§*CLlZl *v*or l3) 

CL  1X1 )  a  0*0 

CLlXX)  a  -X7RCRT-rxrS)-ZC0*fXrSl 

GL<xi)  a  -rxroi-zco*cL(iii*v*oR*r3i-Tv*CLaoi 

CLIXA)  a  -V*CL(XO) 

CLlXX)  a  XNfBT.F  |T)*ZCQ*r(1)-ZCi*CLU1> 

CLCtA)  a  rro)-zco*CL<tO)*iC9*rr(io)-v*or{S))*v*or(4) 

CLIXT)  a  OISRLMH 
40  RETURN 
rNO 


RU*PnUTlNE  EXCITE 

COMMON  /  CONOA  /  RIfOAMMA, ORAV.RO 

common  /  PASO  A  /  8Xl*DXSXl»TMASS.  VNERT.HTaR  (21  > .  ARr  A  (2l )  • 

SZCOE (21)  .OROrT  ill)  «ZlARl21 )  *xt(Zlt»xxSO(|ll* 

X  OVEIOH 12 1 1 .OHASSItl I .ZvT Itl) *ORL 121 1 *ZC9*XNERT. 

X  XZPCRT.OM.MlNKRl.MAXKRI.lNCRCS.ROLnRF 

common  /  T01P  /  vt.MCN.ANSltl.10NrL.riMO.nl 

COMMON  /  MIMO  /  lA.NS.Dri.V.VANOtOMCOA*wAVCN.CM.OtX(21.S)«EACiMA 
COMMON  /  COMO  /  CV(12) *CL(27) .2V.NV.YV.NV.KV 
COMPLEX  2V.Ma.YV.NV.KV 

COMMON  /  BMDa  /  CxrsTl21i.Cxril21)»CXNR<2U.4INMr»1.3ltSRNPrSl.3i 

complex  cxrsr.cxrt.cxNR 

COMMON  /  PP'JOPAM  /  $TORAOEI*ST).YI21).V(21> 

TT  a  U 
M  a  MS 
nx  a  OXT 
MM  a  VAVEN 
CVAM  a  COS(VAnO) 

Span  «  $1N(van9> 

C-01  CALCULATE  VAYE  EXCITATION  AT  tACH  ETaTION 
DO  TC  lal.NS 
**CM  •  -MN*H  (I)*CVAN 
fl<  a  COS(XRCV) 

S«r  a  SIN1XKCV1 

r*v  a  -EXR<-VN*n*AFT 'I)*SECOE It ) )*«A 
XA  a  BSTAR ( I ) *Vn*S VAN/2.0 
ir  e  XA,f 0.0,0  i  00  TO  it 
r*Y  a  tlY*SIN(XAl/Xf 
1?  ir  l  KL.0T.2  )  00  TO  10 

C-02  rolM  VERTICIL  FORCE  COMRQNtNTJ 

ckl  a  QAMMA*0ST£Rrl)-VN«9RAV*ANS:i«n 
xv l  •  vN*c**rAN4(i.2)*EAc-v*o:xrl.ln 

rx  •  I  EKi.*SXK*SKL*CXK)*EXT 

XX  ■  <  r*L»CXK-S*L*SXK»*fIY 

CxrxTlU  a  CMRL* (CX.SX) 

Tr  (  VU.lT.3  )  <10  TO  30 

C-01  ro»M  LATERAL  FORCE  COMPONENTS 

10  XVL  a  OPAV*(RO*APEA(l).ANS(T»3)-VN*ANS:i.5)) 

XKL  a  CM*rANS(I.A)-V*OtX(t.2).VN*V*OIX(1.3n 

r.Y  a  aN*EXY#SVAN 

Cf  a  (  r«L*CXA*SKL#SXN)*EXY 

XX  a  r-EVL*SXV.SrL»CXK)*tXY 

CxCLtI)  ■  CM>>L*rcx.Sx) 

C-0*  FORM  POLL  MOMENT  COMPONENTS 

Ml  a  OPAV*IRO*r9STARrl)**3/12.  -area  r 1 ) *ZI AP ( I ) ) -AMS  1 1 «S) ) 

x  -7C0*rvu 

x«l  •  CM*r-ANau.A. •v*oixrt,3))-7Co*SKL 
ct  •  \  Ml*CXH*S*L*SXK)*tXT 


•  M.tftl  a  CAi'S <PA)*y7*2f STT1 

■*1101  a  ATlN2(PEALlPA).AlM40epA))*S7.2«ST7« 

20  *rTU»N 
END 


SX  m  (-•KL*SX«*S*l*e*'0*C*Y 

rx«*H)  •  :«^*(C«»IU 
Jo  continue 

if  I  «L,«t.7  |  «0  TO  40 

C-OS  INTEGRATE  VtftTICAL  PORCC  ANO  »1TCN  MOMENT 
00  V  !«1«N5 
yd)  •  •CAk.texrjTttn 
it  «m  •  *iM*A(C*r»ran 
C*  •  SJNTIIT.m.y.O*) 
u  »  siN~<ir*M,M,oi) 

7M  ■  f-t»La!CX**ll 
00  M  I-l.Nl 
vdt  •  v Cl >•*! < f i 

31  w<t>  ■  V(t)*xi(t> 

C»  4  -SINT(IT#P,Y.D*) 

4X  •  -IINT  UT«M.M,0*1 

mh  •  cm»l*icx*si> 

|P  (  BJ.LT.3  )  to  TO  SO 

C-CA  tUTCOHAlf  L*T£RAL  PORCEt  YAfc  MOMENT  ANO  *  ' «.L  MOMfNT 

40  no  4X  !•!**•$ 

Yl  I)  •  MALlCAFl.  <!M 
4 1*  Ml!)  •  AlMAfllCiri  (in 
C*  •  S!*T(IT#m,y,OX) 

4x  •  si*-tnT.**,»»,oj) 
v»*  •  cmplx<C*«Sx> 
no  4i  jslihS 
Till  •  YllltXWI) 

4)  Will  •  W(X)«M  d) 

o*  .  sintiit.m.y,o*> 

SX  •  SINTUT.m.w.OX) 

NM  •  C**AL*  <C*»Sl> 

no  4*  !•!«*.$ 

veil  •  »£ALfCA“P  < 1 >  > 

44  Will  •  AlMAft(CX*ft  (t)} 

r*  »  sxnt(it#m.y.o*> 

4X  m  . . . 

mm  •  CMPLX(CXtSX) 

so  return 

fwo 


Su*»rutlNE  WOTlON 

rO****ON  /  TDI*  /  mE.mCN,AnS<71  .101  .KL.KU.IO.XM 
ro  <*pn  /  COwO  /  Cv < 17) iCl  (27)  .Zm.mm.yminm.km 
r  ,«W0N  /  NOT I  /  7A.TA.SA.YA.RA 
OMPlCX  ZA.TA.SA.YA.RA 

COMMON  /  MOTN  t  7HI?l>#7Rl9l).TMI51)*TP<91>.SN(91|,<R(9l).VM|51t. 

X  yp(51)  .AMisjj  ,om> 

CO**lfA  *»«.•. S.T.UtV.M.X.ZM.MM.YM.NN.KHdEN 
mES  • 
w»  •  -yf 

If  I  XL .ST. 7  I  *0  TO  10 
C-01  VP»MCU  MOTIONS  COMMUTATIONS 

•  •  cnmlmcvi  3i«»ts*cv(  |j.mi*cvi  In 
o  •  .CM’ixccvi  Ai-wes«cvi  4).vx*ev<  sn 

•  •  .CMfi*«cvnii-4Es*cv(io»»wx4cvinn 

%  m  C»’lttCv(  4)-*ES»CVl  7>.V**CVI  •  )) 

0/N  4  M#S«4»0 

7 a  •  i2w*s-Mfc*oi/oeN 
TA  •  |M«M«-4»2w)/0CN 
7MIJA)  •  CAASIIA)  *'  * 

7PI10I  #  ATAn2(4EAL(7A>.AIMaS(?A>>M7,HS77* 

Tmijni  •  CAMSITA)*S7.2«S7T0 

»*HA>  •  ATAN2l*CAL(TA)*AlNA9|TA|)*97.2«977« 

*c  I  mu.iT.i  I  no  T0  70 

C*fr7  iATfoai.  -0*10*9  COMMUTATIONS 

in  >  •  cmmlai'u  T>.yts«Cw(  1).m**cl«  in 
o  •  cmpl»<cli  4»-ytS'Cti  4 1 » m**cl (  5n 

•  fc  CMt-CnCLi  *I-VCS*CLI  *>.M**CL(  •»» 

4  *  CKi»n<riti7i-yf;*cui9>.Mx»CLiU>J 

•  •  C**Mt x >CL  H9!*vtS*CL  (13) . « »*Cl  (14|| 

••  .  c»»i.x<CL(i">-»fs«:t.<iM'Vi*ci<!7n 

y  c  CMI*l»ICii7C.4CS*CL<19l(«i*Cwlt0n 
W  •  C«C(.X(ri.t7A>-*CS*CL(2?>tW>«CL<23>! 

»  ■  C“PL*  ICLdM-nFS-fl.  l?9i  ,h*«CL  (»*>  > 

nr*  •  P«  T»  X.  S*  U»  V*  »•  J«  «•  V*  T •  P.  V*  U4  C.  It  s*  0 

4 a  >  |yw*  t •  x*  o*  u»4N»  *-*»•  T«  •-  w*  u*r«.  mn*a  0)/0Cn 

V*  »  1  M*Nfc4  l*y»»  u»  V*  **4N«  U«  Pm  !•  fAYNI/OtN 

•  A  •  c  P*  T»AK.  V*rw«  S*  «•  V  T»TN-  N«NNA  P-KU*  *•  4)/0CN 

SNIIM  •  CAMSiSA) 

4Plt*>  •  ATisJ IACal (SA) • A Imao |SA) )*S7«2SS77S 
vwilft)  •  CAAS(TA)*S7.?SS77« 

vPItA)  •  ATan2{4FA(,  |TA|  *AINA4(tA)  ) *9 7*299779 
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SUBMPUTtwC  HNOSM 

rO"MAN  /  CCNOA  /  PI»«AKWt(MAV»P0 

rpMMON  /  NHQT  /  M0All4)tDTAi0TH«I9iICtlD*ir«Xf >I«»lMtlI*IJ«tTS(9) 
COMMAN  /  HASPa  /  •ML»0IJML«TNAJS,VNCAT,|STA»(21|,A»fA(*n* 
x  sccoeuntoftAFT'in  »z«Ain>n  tiiiin.xiioiii)  • 

X  0«CX«NI|1)  •0HAIK21)  .IdTlin  .MLI2)  ).ZC«.xNfAT, 

X  l2Mf*T.#M.wI^»l,MAXA»I,lNC»CS.»0LnP» 

nOMMAN  /  T0!»  /  wrt«CNtANSllltlO)»KLiRU»10»I« 

COMMON  /  MlMD  /  I A  »N9»0x 1 1 V*  MAN#aOM£9A • pAvCNtCV* 011 (21 iS) tFACtVA 
OOMmAN  /  6*0A  /  CirjT(n)  tClT(.(71)»CIM"(2I),S«MM(9l,3)  tSINPIflt)) 
poMiar*  fxFST,cifL»cxN« 

COMMON  /  MOTI  /  24*TP.S»*yP«»N 

complex  Z»»tp«sp»tp*pp 

COMMON  /  PROOPAN  /  CLlNlti*3) >ec3H(91*2) .SPACE (I0> .SHMIJ) ,fMP(J> , 
l  RMm i  j j  tPHRO)  (HE 1 1 ZPO » Z*CO •  7P0« TP0O *990* 9*00 tVPOtYPOOt MO t 

X  RP0n«CTFST<21*)>(A(4 

r amflE*  «EI(2M0(2P00(TP0(TP00($P0(SP0O(VPO(rPO0«PPO,PP0O( 

X  CTASTtA.B 

C-OI  4fT  lip  CALCULATION  PAPAHCTEPS 
Hfl  •  CnpL*»0#0.-hE> 

JL  •  <<l*9>/4 
,IU  •  IX'J(9)/4 
**m  •  o*i/i. e 
NT  ■  NS-1 

C-07  C ALCVLATE  TOTAL  LOCAL  LOAOlNtS  AT  EACH  STATION 
IP  I  M  .0T.7  I  «0  TO  II 

C-03  vertical  po»ce  component* 

7»D  4  2P4HII 
TIO  •  T4*«fl 

7400  ■  7P0**EI 
ttpn  •  T»o*yti 

no  in  i»i(Ns 

10  CTE4T<t,i)  •  -(AMAS5(X)*ANS(I.in4(2PD0-XIlI)4T»00j  .ANSll.l)* 

*  2»0*V*TPO  •0AAMA*MYAP(I)*<2P-xIII)4T»)  -IanSII.I)* 

X  EAC-V*CI*<I*ln#C?P0-*HI)MT»0»V4TP)  .CXE9TCI) 

IP  I  *U.LT.3  )  SO  TO  1« 

C-OA  LATPPAL  PO»CE  AmO  TORSIONAL  nOm£NT  COMPONENTS 
17  440  m  SP4MEI 
YPO  »  TP4HEI 
DPO  •  PB4HEI 
4P00  •  seo»MEI 
yppn  •  y«o*vEi 

P»OC  •  »PD*Mtl 
AO  14  I41.NS 

ctp4t 1 1 ( 7>  •  •  AMASsii»M:4poo*xim4Y*oo»2Hfm«R»oAi 
x  -ANS(1i3>*tSROO»XI<II#VftOO*2.C4V»YPD>  -IANSCt(41-Y» 

X  DTx (Id) t a (SRQ»Xl <I)4TR0-V*rR*2CS#*R0)  •(AnS<!«9)* 

X  2CO*ANS(t(3l  1*RR0I>  •  ( ANS ( 1 , 10 ) -V*OIX ( 1 ,9)  I «*R0 

X  •C*PL(I> 

CTP4TII *3)  -  -OMASS<I)*(ORLIt)**2**ROD-;«T(t)4(SROO*Xia)«YROO 
1  -0BAVAR»))-04«MA*l4STARa)*  3/17.-A4E4 1 1  >  ♦«7»A*C  I )  « 

>  TC0n«*»  -RROO* t*Nf(  XtT)  •ZC0»(ANS(I»91  «AwS(  I*V>  «7Ct*ANS(  I *3) ) ) 

S  •2C«*CTPST(Id)4lJ  •RRO*(V*lOmi(4l*ZCt«lOIX(t«3l 

4  *01?  < I «9) ♦ZCS’OIX  1 1 #7) ) )  *ZCS* (ANSI I *10) «ANS( I « A) 

9  •reO*ANS(tt4))*ANS(I«S)**OLORP/l»L)*(S*OO*XI(I)4y*O0d.O*V*YRO) 
A  « (ANS(I»9)«ZC»*ANS(I.3n  •(S*O*Xt(I)*YR0-V«vR)« 

7  (ANSir,4>*4-Cfl4ANS(:.4)-V4lOIX(1.4)«ZCS*OIXC«?)n*CXMt|I) 

14  rnNtjNur 

C-C4  (OOP  OVER  STATIl'NS  EOR  SENOINO  MOMENT  CALCS. 

14  WRIT  •  MINKRI 

IP  I  KRIT.OT.O  )  00  TO  IS 

M  #  lt(l)«MH*fl.)-XA) 

AO  tA  14 

1R  XX  4  XI(K*IT).Mh4i1.0.IA) 

C-OS  LOOP  OVER  NUMit*  OP  TvPES  OP  LCAJINQS 

14  no  TA  F4JKJU 

•  •  (C.  0.0.0) 
p  •  (0.0. 0.0) 

IE  (  KRJT.EO.y  )  30  TO  22 

A  •  cTE<T(1«K)/(I*IA) 
p  •  *•(»! (J)-XX) 

TE  (  4RJT.C0.1  I  10  TO  27 
PA  7A  UftPRlT 
A  4  A»CTESttX*K) 

7«  P  •  ••CTESTII.H)M(XT (I)-XX' 


I 


1 

I 

i 


l 


21  *«  ■ 

|R  I  HK.OT.M  )  ao  TO  24 
a  •  a-ctrst(ns,p>/(1.xa> 

4  •  ••CTrSTlM*,*|/()*!*J«(lf <N*)-XX) 

If  t  pk.ot.nt  i  ao  to  24 

00  2*  tRKK.NT 

*  •  <-CTf<Trr«*» 

24  ■  •  o-CTfSTM«Kj*«xi(n-y*s 

24  :r  <  n.ro.3  )  a  •  a 

SMM»P>  ■  CAAJ  (  A  1  *KH 
PMM(«j  .  C*4S(4>»nM 

<H»nn  •  ATANciRCAL  I*)  ,AIMA0(A)1  •110,0/H 

*  m4 (K )  ■  ATAN2(rE*L(R>  tAjMAaO)  >*1«0.0/4I 

30  CONTINUE 

If  <  MAKKfil.ea.HjHAB!  I  CO  TO  3} 

PRINT  4ft,  0**Ct,4,«AlT#(«**M»l),**<AU)»l»lO) 

If  I  K'fT.hf .<NS-IM/l  )  to  TO  3* 

C-OT  4*0*7  WSULT 

31  no  32  K-jL.jy 

4»MM|I0«t<)  •  RMM(K) 

32  *MRI!Ot*l  ■  PMP(«> 

34  tr  {  «OIT.GE,mA»k«I  )  00  TO  3* 

»8TT  •  K»|T*XNCI»eS 

no  To  in 

34  t*  i  JH  *LC*  01  ao  TO  40 
2-1)*  CMffir  SMEAR  A  NO  SENDING  MOMENT  CLOSURE 

no  34  4«1«2 
c lxm.io.ki  ■  o.o 

34  ClA^ilO.ifl  •  0.0 
U*“ll0.3>  •  0.0 
no  o'*  prjl.ju 

*  •  (CTESTIl»A)«CTRST(NS*Kn/M*IAI 

A  •  (CTrST<)*ft)*Xl(l)«CTMT<NS.K>*X!<N$l)/n«lA) 
no  A*  T*2*NT 
A  •  A.CTRST(X.K) 

*0  4  ■  P.CTRST  1 1  **)•*!  ID 

ir  i  k ,i t.3  |  an  to  as 
a  «  « 

no  To  ao 

**  CLSmiIO.ki  *  CA4S(i)«DXt/0I$PL 
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